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ABSTRACT 

The Green Bank Telescope Hll Region Discovery Survey has doubled the 
number of known H II regions in the Galactic zone 343 ° < i < 67° with | | < 1 °. 
We detected 603 discrete hydrogen radio recombination line (RRL) components 
at 9 GHz (3 cm) from 448 targets. Our targets were selected based on spatially 
coincident mid-infrared and 20 cm radio continuum emission. Such sources are 
almost invariably H II regions; we detected hydrogen RRL emission from 95% of 
our target sample. The sensitivity of the Green Bank Telescope and the power of 
its spectrometer together made this survey possible. Here we provide a catalog 
of the measured properties of the RRL and continuum emission from the survey 
nebulae. The derived survey completeness limit, 180 mJy at 9 GHz, is sufficient 
to detect all Hll regions ionized by single 0-stars to a distance of 12kpc. These 
recently discovered nebulae share the same distribution on the sky as does the 
previously known census of Galactic H II regions. On average, however, the new 
nebulae have fainter continuum fluxes, smaller continuum angular sizes, fainter 
RRL intensities and smaller RRL line widths. Though small in angular size, 
many of our new nebulae show little spatial correlation with tracers associated 
with extremely young Hll regions, implying that our sample spans a range of 
evolutionary states. We discovered 34 first quadrant negative- velocity Hll re- 
gions, which lie at extreme distances from the Sun and appear to be part of the 
Outer Arm. We found RRL emission from 207 Spitzer GLIMPSE 8.0 /zm "bub- 
ble" sources, 65 of which have been cataloged previously. It thus appears that 
nearly all GLIMPSE bubbles are Hll regions and that ~ 50% of all Galactic Hll 
regions have a bubble morphology at 8.0 /im. 
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Introduction 



H II regions are zones of plasma surrounding massive O- and B-type stars that are kept 
ionized by the stellar far ultraviolet (FUV) radiation. Stars capable of creating H II regions 
have lifetimes of ~ 10 million years. These nebulae are zero age objects compared to the 
age of the Milky Way and are therefore located at sites of recent massive star formation. 
Hll regions are the most luminous objects in the Milky Way at mid-infrared (MIR) to radio 
wavelengths and can be seen across the entire Galactic disk. Because of their luminosity, 
they are the clearest indicators of the active sites of massive star formation. They are the 
archetypical tracers of Galactic spiral structure. Their chemical abundances, moreover, indi- 
cate the present state of the interstellar medium (ISM) and reveal the elemental enrichment 
caused by the nuclear processing of many stellar generations. They provide unique and 
important probes of billions of years of Galactic chemical evolution (GCE). 

Knowing the Galactic distribution, kinematics, and physical properties of star clus- 
ter/HlI region/PDR/molecular cloud complexes informs a host of research areas, including 
the merger history of the Milky Way, the structure of the Galaxy, the physics of star for- 
mation, and the evolution of the ISM. At the most basic level, a census locates new sites of 
massive star formation and paves the way for future studies at different wavelengths. Indi- 
vidual Galactic Hll regions are astrophysically important objects that reveal details of the 
impact of the star formation process on the ISM. 

Hll region RRL spectra yield the velocity of the nebulae, which can then be used 
to derive their kinematic distances . The distances place Hll regions onto a G alactic map 



that traces Galactic structure (e.g. iDownes et al.lllQSOt lAnderson &: Baniall2009[ ). The RRL 



line width measures a combination of the turbulent and thermal energies in the Hll region 
plasma and is a parameter needed to derive the electron temperature of Hll regions. The 
electron temperature is a proxy for nebular metallicity and hence is key to understanding 
the M ilky Way metallic i ty gradient, which provides an irnportan t constraint on models for 
GCE dwink et al.lllQSsi Ishaver et al.lEgSsI : lOuireza et~aDl2006b[ l. 



Radio continuum observations of Hll regions measure the free-free thermal emission 
and give the radio flux and angular size. When combined with a distance estimate, these 
parameters give a luminosity and a physical size. The former constrains the spectral type of 
the ionizing star and the latter is related to the evolutionary stage of the H II region. Finally, 
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with a large volume-limited sample, we may better estimate the number of Hll regions in 
the Galaxy and derive their luminosity function. Studies of external galaxies have found 
that the shape of the Hll region lumi nosity distribution is related to the morphology of the 
host Galaxy (see lOey fc Clarkelll998l and references therein). By deriving the slope of the 
luminosity distribution of Galactic Hll regions, we may get a better understanding of the 
morphology of the Milky Way. 

Modern Galactic Hll region s u rveys began with studies of the Palomar Sky Survey 
optical plates ( 1Sharplesslll953l . 119591 ) . iGumI ( 1l955l ) found 85 nebulae vis ible from the Souther n 
hemisphere. The Sharpless and Gum surveys were soon augmented by lRodgers et al.l ( Il960l ). 
who cataloged 182 optical Hll regions in the Southern sky. (Some of these are also Sharpless 
nebulae.) Extinction within the Galactic plane, however, limits the detection of optical 
nebulae to within a few kpc of the Sun. At visible wavelengths the Galactic census of Hll 
regions thus ends at the local Solar neighborhood. Only a census based on observations at 
longer wavelengths, where the ISM is optically thin, can locate Hll regions on a Galactic 
scale. 

Hll regions emit at radio wavelengths because of their thermal free-free continuum 

(Bremsstrahlung) and radio recombii iation line emission. Radio continuum observations of 

the Milky Way began in the 1 950s (see Westerhoutj 1958a : Piddington 1951; Hanbury Brown &: Hazard 
19531 : IScheuer fc Rylelll953l : iHaddock et al.lll954l ). Optically visible Hll regions were also 
bright at radio wavelengths. Moreover, most of the bright, discrete radio continuum sources 
in the Galactic plane turne d out to be thermall y emitting H II regions. For example, of the 
74 continuum sources in the IWesterhoutI ( 1l958bl ) catalog (not counting W82, the Moon, and 
W24, the Galactic center), 55 (75%) are thermally emitting Hll regions, 11 are non-thermally 
emitting supernova remnants (SNRs), and 7 are non-thermally emitting extra-galactic active 
galactic nuclei (AG N). Of the 55 Westerhout Hll region sourc es, 30 are found in the catalogs 
of optical nebulae (ISharplesslll953l . 1 19591 : iRodgers et al.lll960l ). 



Radio recombination lines from Galactic H II regions were discovered by lHoglund fc Mezger 



( 119651 ) who detected H109a emission from M17 and Orion A. Because the Galactic ISM is 
optically thin at centimeter wavelengths, RRL surveys were able to discover large numbers of 
HlI regions distributed throughout the en tire Galactic disk. These p i oneering surveys wer e 
done in the 1960s, 197 s, and 1980s bv iMezger fc Hodundl fll967 ): Iwilson et al.l /ll97oh : 



Reifenstein et al.l jl97ol ): IPownes et all Jl98ol l: ICaswell fc HavnesI Jl987h : Ihockmanl jl989 



the latter hereafter referred to as L89). These surveys gave important insights into Galactic 
structure and the spatial distribution of massive star formation. Particularly noteworthy 
was the discovery of a metallicity gradient across the Galactic disk, made apparent by RRL 



measurements of H II region electron temperatures ( iWink et al.l Il983l : IShaver et al. 



1983 
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Quireza et al.ll2006b( ). By the time of the L89 survey, however, almost all of the reasonably 
strong radio continuum sources had been observed; the iLockman et al.l ( 119961 ) study of "dif- 
fuse" H II regions was the last large angular scale survey for discrete H II regions made using 
cm-wavelength RRLs as tracers. 

Most recent work on Galactic H II regions has focused on the ultra-c ompact (UC) class of 
nebulae, whi ch nominally are in an earl y phase of H II region evolution. IWood fc Churchwell 
( 1l989a[ ) and iHughes fc MacLeodI (119891 ) used the colors (fl ux ratios) of sources f rom the In- 
frared Astronomical Satellite [IRAS] point source catalog ( jBeichman et al.lll988l ) to identify 
candidate UC Hll regions. Hll regions are bright at infrared (IR) wavelengths largely be- 
cause radiation from the central star (or stars) is absorbed by lo cal dust, which thernaally re- 
emits at IR wavel engths. Later surveys of the radio continuum (IWood fc Churchwell! Il989b 



Kurtz et al.l Il994l ) and RRL emission (lAraya et al.l l2002l : IWatson et al 



2003 



Sewilo et al. 



20041 ) of these targets showed that IRAS color criteria can be used to identify UC Hll 



regions. 

Despite these efforts, the census of Galactic Hll regions was clearly incomplete. The 
advent of modern high-resolution. Galactic-scale infrared and radio surveys (see §2]), coupled 
with the unprecedented spectral sensitivity of the NRAO Green B ank Telescope (GE T), has 
made the GET Hll Region Discovery Survey (HRDS) possible ( lEania et al.l l2010l ). The 
new radio a nd IR surveys have greatl y increased our ability to identify Galactic Hll region 
candidates (iGiveon et al.l l2005l . |2008| ). The increases in sensitivity due to advances in in- 
strumentation and telescope construction allow us to detect Hll region candidates in radio 
continuum and RRL emission. In particular, the GET's unblocked aperture, active surface, 
and Autocorrelation Spectrometer give it unprecedented spectral sensitivity at centimeter- 
wavelengths. 

The GET HRDS extends over 168 square degrees, covering the zone from 343 ° < i < 
67° with \b \ < 1°. It dete cted 603 discrete hydrogen RRL components at 9 GHz (3 cm; 
"X-band") from 448 targets (lEania et al.ll2010l ). Here we provide a catalog of the measured 
properties of the RRL and continuum emission from the GET HRDS nebulae, compare these 
properties to the previously known census of Galactic Hll regions, and identify astrophysi- 
cally important classes of sources discovered by the HRDS. These include populations of Hll 
regions that lie at extreme distances from the Sun in the Outer Arm and that are associated 
with Spitzer GLIMPSE 8.0 /xm "bubble" sources. 



- 5 - 



2. The HRDS Target Sample 



We assembled our target list from the following multi-frequency, large solid angle Galac- 
tic surveys: th e NRAO VLA Galactic Plane Survey at 21cm Hi and conti nuum (VGPS: 



Stil et al.M2006h . the NRAO VLA Sky Survey at 20 cm continuum fNVSS: Condon et al. 
IQQsh . the Southern Galactic Plane Survey at 21 cm H I and continuum (SGPS: Haverkorn et al 



200d ]. the VLA M AGPIS at 20 cm c ontinuum Jflelfand et al-lbood ). and the Spitzer 2Aij.m 
MIPSGAL survey JCarev et allboooh . 



The HRDS candidate Hll regions are targets that: (1) have spatially coincident 24 /xm 
MIR and 20 cm radio continuum emission of a similar morphology and angular extent, (2) 
have an extrapolated 9 GHz flux density of at least 70 mJy, (3) are not known to be an H II 
region because of a previous RRL detection, and (4) are not known to be a planetary neb- 
ula (PN), supernova remnant (SNR), pulsar wind nebula (PWN), Active Galactic Nucleus 
(AGN), or a luminous blue variable star (LBV) . Sources having both MIR and radio con- 



tinuu m emission are likely emitting thermally (IHaslam fc Osbornd llQSTt iBroadbent et al. 



19891 ). and are therefore probably either Hll regions or planetary nebulae. The plasma ion- 
ized by the FUV radiation from the exciting star(s) gives rise to free- free thermal emission at 
cm-wavelengths. Warm dust in the nebula absorbs the stellar radiation as well and re-emits 
in the MIR. Furthermore, small dust grains, which are stochastically heated by the stellar 
flux, also emit in the MIR. 

For the radio continuum we primarily use data from the VGPS. This survey measured 
21cm Hi line emission, but produced 21cm continuum maps from spectral channels with 
no line emission. The VGPS extends from 17°5 < £ < 67 ?5 with latitude coverage varying 
between 1?3 and 2?6 in \b \ and it has an angular resolution of 1'. The VGPS combines 
data from the GBT to fill in the zero-spacing information and is therefore sensitive to both 
small and large-scale emission features. For the longitude range 358° < ^ < 17?5, we use the 
NRAO NVSS 20 cm continuum data. The NVSS has an angular resolution of 45" and covers 
the entire sky north of —40° Declination. It does not include zero-spacing information and 
therefore some larger, diffuse emission regions are not detected. For 343° < i < 358°, we use 
the continuum data from the SGPS, which is the southern counterpart to the VGPS. The data 
were obtained in a similar manner using ATCA for the 21cm interferometry measurements 
and the Parkes telescope for the zero-spacing information. The resolution of the SGPS is 
100" and it extends from 253° < i < 358° with \ b \ < 1°. Finally, to get the best positions 
for targets in regions of complex continuum emission, we use the MAGPIS 20 cm data. 
MAGPIS was made with the VLA and extends from 5° < i < 48° with | 6 | < 0?8 at an 
angular resolution of ~ 5". Although MAGPIS has the zero-spacing information from the 
Effelsberg 100 m telescope, visual inspection shows that diffuse emission detected with the 
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VGPS is often not seen by MAGPIS. We prefer to use the VGPS, NVSS, and SOPS since 
their resolutions are more comparable to that of the GBT at X-band, and only use MAGPIS 
to disentangle complex emission regions. 

For the MIR emission we use MIPSGAL 24 n m data, which w ere obtained with the 
Multihand Imaging Photometer for Spitzer (MIPS: iRieke et al.ll2004l ) on the Spitzer Space 
Telescope. MIPSGAL covers 293° < £ < 67° with | 6 | < 1° at a resolution of 6" at 24 /im. 
Due to the low 24 /im optical depth of the intervening ISM, MIPSGAL can detect Hll 
regions on the far side of the Galaxy. Hll regions are bright in the 24 /im MIPSGAL band 
for two reasons. First, there is thermal emission from dust grains spatially coincident with 
the GLIMPSE 8.0 /xm emission (see below) produced by the Hll region photo-dissociation 
region (PDR). There is also emission that is spatially coincident with the ionized gas. This 
emission is most likely produced by small dust grains that are stochastically heated by 
absorbing ultraviolet photons from the exciting star(s). The grains have a temperature of 
> 100 K and they are not in thermal equilibrium. At 24 /im the flux of both compo nents is 
roughly equal for Hll regions surrounded by 8.0 /im bubbles ( iDeharveng et al.ll2010l ). 



Our analysis here also uses 8.0 /xm data from the Ga lactic Legacy Infrared Mid-Plane 
Survey Extraordinaire f GLIM PSE: iBeniamin et al.M2003l ). which were obtained with the 
Infrared Array Camera (IRAC: iFazio et al.ll2004( ) on the Spitzer Space Telescope. GLIMPSE 
has the same Galactic coverage as MIPSGAL, but with a resolution of ~ 2" at 8.0 /im. In 
addition to the re-radiated emission from heated dust grains, the IRAC 8.0 /im band contains 
emission from polycyclic aromatic hydrocarbons (PAHs). These molecules fluoresce when 
they absorb ultraviolet radiation, and are thus an excellent tracer of H II region PDRs. 

Criterion T. Spatially Coincident 20 cm and 24fim Emission To find Hll region candidate 
targets we seek sources with matching 20 cm and 24 /xm emission. The MIR flux associ- 
ated with Hll regions is or ders of magnitude greater than the radio continuu m flux (e.g. 



Wood fc Churchwelllll989bf ). Since the MIPSGAL 5 a sensitivity is 1.7mJy fICarev et al. 
20081 ). we expect that all Hll regions with 9 GHz fluxes > 70mJy (see Criterion 2 below) 
should be easily detected by MIPSGAL. 

To find our H II region candidates we visually compare the 20 cm radio continuum emis- 
sion at a particular target location with the 24 /xm emission at the corresponding position, 
searching for sources with spatially coincident radio and MIR emission that also have a simi- 
lar angular extent. We did not automate our method because of the wide range of H II region 
morphologies, together with the complex environments in which H II regions are found. Au- 
tomated detection methods would be extremely difficult to implement, and would likely fail 
in complex regions. Visual inspection ensures that the target catalog includes both small, 
compact sources as well as larger, diffuse regions that automated detection might confuse 
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with the background. To search for Hll regions, we use the DS9 softwar^ ( Jove fc Mandel 



20031 ) to ahgn the radio and MIR images so that their astrometry matches. We then repeat- 
edly bhnk between these images in order to find the nominal position of the target, which is 
the visually determined peak of the radio emission. In order to find all possible targets, we 
made this analysis for each field at least three times. Each pass through the data added new 
targets. By the third pass only a small number were added. These were mostly faint sources, 
but some came from complex fields. After the list was complete, we made a final pass to 
refine the nominal target positions. In addition to this "blind" s earch, we also inspecte d by 



blinking the p ositions of sources fouii d in na olecular CS (2 — t- 1) flBronfman et al.lll996[ ) and 



IR "bubble" (IChurchwell et al.l l2006l . 120071 ) catalogs. These objects are prime Hll region 
candidates. 

Criterion 2: Flux Density The candidate targets must have an extrapolated 9 GHz flux 
density greater than 70 mJy. This is the flux of the thermal continuum radiation from H II 
regions ionized by a single 0-9 star at a distance of 20kpc (see §5.6.ip . There is no extant 
continuum survey at 9 GHz with the sky coverage of the HRDS so we must get the fiuxes of 
target H II regions by extrapolating source fiuxes from radio continuum surveys made at other 
wavelengths. We assume that the sources are optically thin at 20 cm wavelength, so Sy oc 
^y-o i. We use our Kang softwar^ to estimate the 9 GHz continuum fiux of each candidate. 
Our Kang photometry places a circular aperture on the source and an annulus around that 
zone for the background. To be detectable in our survey, optically thin 20 cm H II region 
candidates must then have VGPS or NVSS flux densities > 70 mJy x (1.4 GHz/9 GHz)-°-i ^ 
85mJy. We delete all sources not meeting this flux threshold from the target list. 

Criteria 3 & 4- Removing Previously Known HII Regions and Other Objects Our Hll region 
candidate criteria can be met by a number of astrophysical objects other than Hll regions. 
Active galactic nuclei (AGN), supernova remnants (SNRs), pulsar wind nebulae (PWNe), 
luminous blue variables (LBVs), and planetary nebulae (PNe) can all contaminate the target 
sample. Using the SIMBAD databas^ we remove such contaminating sources by searching 
for previously known objects that are located within 5 ' of the nominal position of all the 
HRDS targets. By correlating the candidate sources with the SIMBAD database in this way, 
we also remove known H II regions with measured RRL emission. 

It is certainly true that the radio continuum can result from a mixture of free-free 
(thermal) and synchrotron (non-thermal) emission. Non-thermal emitters such as AGNs, 



^ ht tp : / /hea- www . harvard . edu/ RD /ds9/ 
^http://www. bu.edu/iar/kang/ 
http://sinibad.u-strasbg.fr/sinibad/. 
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SNRs, and PWNe, however, cannot contribute many contaminating sources to the sample. 
Fuerst et aLl ( 119871 ) showed that one can discriminate thermally from non-thermally emitting 
objects by using the ratio of the infrared to radio fluxes. The IR/radio flux ratio for Hll 
regions is typically ~ 100 times larger than that for non-thermally emitting SNRs, so it is 
easy to differentiate between the two by eye. Our visual inspection of the MIR and radio 
images should eliminate all the SNRs from our target sample. Furthermore, since non- 
thermal emitting sources do not emit RRLs, even if they were on the target list they would 
not appear in the HRDS catalog of detections (see §4.ip . Sources showing coincident MIR 
and cm-wave continuum emission almost invariably are thermally emitting: 95% of our 
sample targets show hydrogen RRL emission with line to continuum ratios of ~ 10^^ which 
together suggests that our targets are emitting thermally. 

Since PNe emit thermally, they are by far the largest source of contamination. Cor- 
relation with SIMBAD helps remove previously known PNe, but there might be many un- 
cataloged PNe that do not appear in SIMB AD. PNe do, however, have broader lines than 
Hll regions because of their expansion (see iGaray et al.lll989l : iBalser et al.l 119971 ). which in 
principle allows the two populations to be separated ex post facto. Based on their line widths 
and compact IR emission, a small number of the HRDS sources are probably uncataloged 
PNe (see M) . 

Because of the emission from their ejected material, LBV and other evolve d stars with 



ejected nebulae may a l so appear bright at MIR and radio wavelengths (see IStephenson 



1992; Clark et al. 2003; Mizuno et al 



2010l for their MIR emission characteristics). Most 



LBVs are removed by correlation with the SIMBAD database using the nominal HRDS 
target positions. We also identify, and then remove, LBVs by visual examination, since they 
appear as nearly circular objects in MIPSGAL with very thin emission rings, whereas Hll 
regions generally show less symmetry and broader PDRs. 

In the HRDS survey zone we find over 1,000 Hll region candidates that have spatially 
coincident MIR and radio continuum emission with a similar morphology and angular extent, 
criterion (1). The flux cutoff of 70 mJy at 9 GHz, criterion (2), reduces this number to ~ 600 
sources, and removing known Hll regions and contaminating sources, criteria (3) and (4), 
gives a final H II region candidate list of 470 targets. Images of the MIR and radio continuum 
emission for four typical HRDS targets are shown in Figure [H These images are 5' on a side 
and show GLIMPSE 8.0 /im (left column), MIPSGAL 24 /im (middle column), and VGPS 
21cm (right column) emission for each target. Evident is the similar angular extent of these 
targets. The environment surrounding the targets can be quite complex. The "bubble" 
morphology seen in the GLIMPSE images is a common feature of our targets (see §4.3p . 
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3. Observations and Data Reduction 

Our observations were made with tlie GBT 100 m telescope from June 2008 tlirougli 
October 2010. For each candidate Hll region, we interleaved spectral line (using the Au- 
tocorrelation Spectrometer [ACS]) and continuum (using the Digital Continuum Receiver 
[DCR]) observations. Interleaving observations in this way gives a RRL line-to-continuum 
ratio that is nearly independent of calibration uncertainties and weather conditions. We fo- 
cused the telescope and established local pointing corrections every few hours using standard 
X-band pointing sources. 

The calibration of the intensity scale was measured to be within 10% for the RRL and 
continuum data. The intensity scale was determined by injecting noise of a known magnitude 
into the signal path at a rate of 1 Hz and 10 Hz for the RRL and continuum observations, 
respectively, with a duty cycle of 50%. The antenna temperature for the line and continuum 
data was set using the measured calibration intensity, Tcai, of the noise source. The value of 
Teal varies with frequency and was measured in the lab at 50 MHz intervals using hot and cold 
loads. We confirmed that the Tcai values are accurate with an uncertainty of less than 10% 
by periodically observing the flux density cahbration sources 3C 147 and 3C 286 under good 



weath er conditions near transit. The flux densities fo r these calibrators are from lPeng et al. 



(I2OOOI ) and we assume a telescope gain of 2 K Jy^^ (iGhigo et al.ll200ll ). We did not make 
any corrections for atmospheric opacity or telescope gain as a function of elevation. Under 
most weather conditions opacity corrections are typically a few percent and at 10 GHz any 



elevation gain corrections are less than 5% (iGhigo et al.ll200ll ). Overall, the intensity scale 



in each ACS band is consistent at the ~ 12% level. 

To maximize the power and flexibility of our analysis we wrote a large suite of IDL pro- 
cedures rather than use any of the standard single-dish radio astronomy software packages. 
We use our own TMBIDL software package to analyze both the line and continuum GBT 
data. Written in IDL, TMBIDL inspired NRAO to develop their GBTIDL softwarffl, which 
is now the main data analysis tool for the GBT. We find, however, that TMBIDL is much 
more fiexible for our needs. We wrote extensive IDL code to analyze ACS and DCR data 
both in real time and also for post data acquisition. This HRDS software led to the latest 
V6.1 version of the TMBIDL package (Bania 2010, private communication). 



*See http://www.bu.edu/iar/tmbidl 
^See http://gbtidl.sourceforge.net/ 
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3.1. Radio Recombination Lines 

The sensitivity of the GBT and the power of its ACS together made t he HRDS pos - 



sible. To achieve high spectral sensitivity, we used techniques pioneered by iBalserl ( 120061 ). 
who realized that there are 8 Hna RRL transitions, H86a to H93a, that fall within the 
instantaneous bandwidth of the GBT X-band receiver. They can all be measured simultane- 
ously by the GBT with the ACS. The H86a transition, however, is spectrally compromised 
by the H 108/3 line at nearly the same frequency so it cannot be used for the HRDS. We 
thus get 14 independent spectra per OffOn total power observation pair (7 transitions x 2 
orthogonal polarizations each). We identify the strongest spectral feature in each band as 
an Hna transition because there are no other strong atomic, molecular, or unidentified lines 
known in these frequency bands. For sources with strong continuum emission we may also 
detect the helium and carbon RRLs. 

The RRL spectra were gotten by using the ACS in total-power, position switching 
mode with On- and Off-source integrations of six minutes each (hereafter a "pair"). The 
Off-source position observation tracked the same azimuth and zenith angle path as the 
On-source scan. The ACS was configured as 16 independent 50 MHz wide spectral bands. 
Each was 9-level sampled by 4096 channels. The typical velocity resolution of any band is 
~ 0.4 kms~^ per channel. Because the target RRL velocities were unknown, all ACS bands 
were tuned to a center LSR § velocity of 42.5 kms~^ for targets in the C range from 10° 
to 67°. For targets in the I = 343° to 10° range the center velocity was 0.0 kms~^. The 
50 MHz ACS bands span ~ 1, 700 kms~^ at X-band, which greatly exceeds the LSR velocity 
range of objects gravitationally bound to the Milky Way. We can thus, in principle, detect 
RRL emission from any Galactic Hll region. 



Balserl (120061 ) showed that all these transitions can be averaged together to improve 
the RRL signal-to-noise ratio, thus giving an extremely sensitive X-band ( H n a ) composite 
nebular spectrum. Because the energy level spacings between adjacent RRLs for n > 50 are 
similar relative to the ground state, the line parameters — intensity, line width, and velocity 
— are nearly identical. For example, the classical oscillator strengths for the H50a and 
H51a RRLs differ by only 2%, so these lines should have comparable intensities (Menzel 
1968). Thus at the high energy levels of our observations, to first order all transitions have 
the same intensity. 

We can thus average all the 14 independent spectra together and get an extremely 



^The RRL velocities here are in the kineniatic local standard of rest (LSR) frame using the radio definition 
of the Doppler shift. The kinema tic LSR is defined by the solar motion of 20.0 km s~^ toward {a, 6) — 
(18'', +30°) [1900.0] |Gordonlll97fih . 
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sensitive ( H n a ) composite RRL spectrum for a source. Three steps are required to average 
the individual Hna RRLs into a source composite (Hna) spectrum: (1) the velocity scale 
must be re-gridded to a common resolution; (2) the RRL spectra must be shifted to match the 
Doppler tracked H 89 a RRL; and (3) the spectra must be aligned in velocity and averaged 
together. The first step is required because the velocity resolution is a function of sky 
frequency, which varies with the different RRL transition frequencies. We therefore re-grid 
the H 88 a to H 93 a RRL spectra to the velocity scale of the H 87 a spectrum, which has the 
poorest spectral resolution. The second step is necessary since in the GBT system only the 
H 89 a spectral band is properly Doppler tracked. For the third step, we align the velocities 
and the average spectra in the standard way with a weighting factor of tintg/T^yg, where tintg 
is the integration time and Tgys is the total system temperature. 

We smooth the ( H n a ) spectrum with a normalized Gaussian over 5 channels to give 
a spectrum with a velocity resolution of 1.86 kms~^, which compares favorably with the 
~ 25 kms~^ FWHM line widths that are typical for RRLs from Galactic Hll regions. We 
then remove a polynomial baseline, which is usually a second order fit. Finally, we fit 
Gaussians to this final composite (Hna) spectrum, which gives us the LSR velocity, line 
intensity and line width for each emission component. The baseline removal and Gaussian 
fitting procedures are not automated, but rather are done individually for each source. 

This observing technique, coupled with the sensitivity afforded by the GBT's aperture, 
gave us unprecedented spectral sensitivity per unit observing time advantage compared with 
all previous cm-wavelength RRL surveys of Galactic Hll regions. The vast majority of our 
targets took only a single OffOn total power pair to detect RRL emission and estabhsh an 
accurate source LSR velocity. Typically, the r.m.s. noise for observations consisting of a 
single pair is ~ 1 mJy. For some of the weaker sources, and for sources for which we wanted 
to measure the H/He line ratio more accurately, we made additional observations. 

Over the frequency span of the H 87 a to H 93 a RRLs, the GBT HPBW beam size varies 
by 20%. We therefore expect that source structure will make the various transition intensities 
differ by more than what is theoretically expected. We observed this effect for bright, 
compact Hll regions that we used as calibrators, such as W3, where the RRL intensities 
approximately scale with the beam area, as is expected for a point source. For the HRDS 
catalog, we did not scale to a common beam size; the (Hna) intensities are averages over 
the beam size range of 73 " to 89 

We show in Figure [2] ( H n a ) spectra for the Figure [T] H II regions. These are typical 
HRDS (Hna) spectra. The G 31.727+0.698 nebula, however, is one of our weakest detec- 
tions. Both hydrogen and helium RRLs can be seen in the G 31.470— 0.344 spectrum. We 
detect two RRL components along the hne of sight to G 30.378+0.106, one hkely due to W43 



(see 



3.2. Radio Continuum 

We used the DCR to measure the continuum flux for our HRDS targets by slewing the 
telescope across the sky while sampling the total power every 100 ms. We made cross scans 
centered on the nominal source position consisting of four total-power integrations: forward 
and backward in Right Ascension, RA, and forward and backward in Declination, Dec. Each 
scan was 20' in length and we slewed at 20" per second. The DCR was centered at 8665 
MHz with a 320 MHz bandwidth. 

The continuum data can be considerably more complex than the RRL data due to a 
combination of source confusion and sky noise. A typical observation of an HRDS target 
starts with a continuum measurement followed by a spectral line pair. For some sources, 
however, we have many such continuum scans. We thus average all the forward and back- 
ward scans (after flrst flipping the backwards scans). This average usually blends different 
observing epochs and weather conditions so it tends to minimize the effects of sky noise. 
We remove a polynomial baseline, usually a second order flt, from the RA and Dec average 
scans and then flt Gaussians to them. This analysis cannot be automated and was done on 
a source by source basis. 

The inner Galaxy has a high density of sources that are not always well-separated by 
the 82" GET beam. We attempt to resolve this confusion by visually examining VGPS, 
NVSS, or SGPS images. These data helped us to separate the target's emission from that of 
other contaminating sources and sky noise. Each source was analyzed individually and the 
polynomial baseline was removed based on what we judged was the baseline zero intensity 
level of the scan. Some sources showed multiple Gaussian components. Sometimes all 
components can be associated with emission from the HRDS target. In other cases one or 
more Gaussian components are due to other, unassociated radio continuum sources. The 
auxiliary radio continuum data from the various sky surveys provided vital extra information 
that we used to decide which components are due to the HRDS target. 

Clearly the resolution of source confusion is a subjective process. Despite our care, the 
flts to the radio continuum data are, especially for weak sources, uncertain and in some cases 
no flt is possible. Figure [31 which shows the continuum scans for the four Figured] nebulae, 
illustrates the complex nature of these data. The thick lines show the Gaussian flt to the 
HRDS component of each scan. For sources with simple emission structure which is also 
well-centered in the GET beam, we expect the fltted heights in the RA and Dec scans to be 
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of equal intensity. Additionally, all sources should have fitted sizes at least as large as the 
GBT beam. When these conditions are not met, it can be a sign of an inaccurate position 
or poor data quality. Continuum parameters for such sources should be used with extreme 
caution. 

For each source we analyze the RA and Dec continuum emission scans independently. 
For targets with simple structure we fit a single Gaussian component, which gives a peak 
intensity (in mK) and full width at half maximum (FWHM) angular size (in arcsec) for each 
RA and Dec scan. For many of the HRDS sources, however, the continuum emission is too 
complex to be modeled by a single Gaussian. For these cases, we fit the continuum scan 
with multiple Gaussian components, each with its own peak intensity and FWHM angular 
size. The sum of these components gives an estimate for the total source integrated intensity 
(units are [mK arcsec]). We then crudely approximate this complex emission with a single 
Gaussian model whose peak intensity and FWHM angular size together produce the same 
integrated intensity as the sum of the multiple components. The FWHM angular size of this 
fictitious single component is set by the maximum separation of the FWHM position extrema 
of the multiple components with respect to the nominal target position. This angular size, 
together with the integrated intensity determined by the multiple components, defines the 
peak intensity of the fictitious single component. 

We use the Gaussian components to estimate the source integrated fiux density at 
8665 MHz. We take the arithmetic average of the RA and Dec intensities to define an 
average peak intensity and use the 2K per Jy GBT X-band gain to convert this into the 
peak fiux. Assuming the source continuum emission is a two-dimensional Gaussian in RA 
and Dec, the integrated source fiux. Si, is: 

^RA \ / ^Dcc 



where Sp is the peak fiux, 9b is the GBT beam size at X-band, and ^ra, ^per are the FWHM 



angular sizes derived from the RA and Dec scans (see lKuchar fc Clarklll997l ). 



4. The H II Region Discovery Survey Catalog 

We detect radio recombination line and continuum emission from 448 of our 470 targets, 
over 95% of our sample of H II region candidates. The HRDS targets with detected hydrogen 
RRLs are listed in Table [1], which gives their Galactic and Equatorial ( J2000) coordinates. 
For brevity, we provide here only data for the first 10 sources in the HRDS catalog. These 
stubs show the information given in each catalog table; the complete HRDS catalog data are 
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provided in the electronic version of this paper. Table [J also provides information about a 



source's a ssociation with metha nol masers (e.g., iPestalozzi et al.ll2005l ) and with CS 2 — )■ 1 



emission (iBronfman et al.lll996l ). together with a characterization of the source morphology 
as seen in Spitzer GLIMPSE 8.0 /xm images. Methanol maser and CS emission are both 
produced in dense gas regions. They are associated with the earliest phases of massive 
star formation; these data therefore give some indication of the evolutionary state of HRDS 
targets. The method we use to associate the maser and CS positions with the HRDS targets, 
as well as the implications of these associations for the HRDS targets are discussed in §4.51 
Sources observed but not detected in CS are listed as "ND". The Table [1] maser velocity is 
that of the emission component with the highest intensity. 

Images of some example HRDS targets are shown in Figure HI They are made from 
Spitzer survey data: MIPSGAL 24 /im (red), GLIMPSE 8.0/im (green), and GLIMPSE 
3.6 yum (blue). HRDS targets are often found in complex IR environments and they have 
a wide range of sizes and morphologies. Figure H] shows examples of all the different mor- 
phological classifications found in Tabled! For nearly all the angularly resolved HRDS Hll 
regions, the 24 /zm emission lies interior to the 8.0 /im emission. This appears to be the MIR 
signature of Galactic H II regions. Three-color Spitzer images such as Figure S] can thus be 
used to find H II region candidates with high reliability. The most striking examples of this 
phenomenon are the GLIMPSE "bubbles" (see §4.3p . 

We detected 603 discrete hydrogen RRL components from the 448 HRDS targets. The 
RRL data are given in the Table |2] stub, which lists the source name, the Galactic longitude 
and latitude, the line intensity, the FWHM line width, the LSR velocity, and the r.m.s. 
noise in the spectrum. The errors given in Table |2] for the fitted line parameters are the 
1 a uncertainties from the Gaussian fitting. For sources with multiple velocity components 
detected along the line of sight, the source names are given additional letters, "a", "b", "c", 
or "d" , in order of decreasing peak line intensity. 

A stub of the HRDS radio continuum catalog data is given in Table |3l which lists the 
source name, the Galactic longitude and latitude, the peak intensity in the RA and Dec scan 
directions, the FWHM angular size in the RA and Dec directions, and the integrated flux 
density derived using Equation [T] For sources better approximated by multiple Gaussians, 
the angular size listed is the maximum extent of the multi-component composite source and 
the intensity is that for a single Gaussian component, which is determined by this width 
and the integrated flux density (see § l3.2p . Such "complex" objects are identified by a "C" 
label in the Notes column of Table El The errors given for the continuum parameters are the 
1 a uncertainties in the Gaussian fits, propagated as necessary for derived quantities using 
Equation [H 
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The main goal of the HRDS was to detect new Galactic Hll regions. To maximize the 
number of observed objects we therefore used standard observing procedures throughout 
the survey and did not customize the observing for specific sources. For example, the peak 
continuum emission observed with the GBT was typically offset relative to the nominal 
target position, but we did not refine these coordinates, i.e., we did not peak-up on the 
source relative to the GBT's beam. Moreover, for some targets it was difficult to establish 
the continuum zero-level given the complexity of the continuum emission structure on the 
sky and time variable sky noise. Source confusion or cases where a very weak target is quite 
near to a strong source can compromise the accuracy of the continuum measurement. In 
some sources longer continuum cross scans might have yielded more accurate baselines, but 
our goal was to produce a new catalog of Galactic Hll regions not to provide the highest 
quality continuum data. Much of the continuum data in the HRDS catalog is therefore 
not suitable for deriving very accurate physical properties of these Hll regions, such as the 
electron density and temperature, the excitation, etc. A subset of sources, however, have 
single component continuum peaks that lie within 10 " of the nominal target position. These 
nebulae, which are flagged with a "P" label in the Notes column of Table [3l are our highest 
quality continuum data. They are suitable for deriving physical properties. 

The HRDS catalog, however, is an excellent resource for Galactic structure studies 
since the nebular positions and LSR velocities are measured with high accuracy. Each RRL 
detection yields the nebular source velocity. The nebular velocity, when combined with the 
(£, h) location, gives the source Galactocentric position and kinematic distance from the Sun. 
The continuum measurements of the thermal free-free emission give the nebular total flux, 
angular size, and, crudely, some knowledge of the source morphology. 

4.1. Sources Not Detected 

We did not detect RRL emission from 22 of the candidate H II region targets. These are 
typically < 1 mJy upper limits, which is the sensitivity of an ( H n a ) composite spectrum 
made from a single observing pair. These non-detections are listed in Table HJ which gives the 
source name, the GLIMPSE 8.0 /im morphology (using the Tabled] classifications), comments 
on each source, and our best estimate of the object's type: Hll region, PN, AGN, or PWN. 
Nearly half of these non-detections, 10 of the 22 targets, may be PNe or emission line stars. 
They are unresolved in the IR and the radio continuum surveys we used to compile our 
target hst. 

Five sources, 23% of the sample of non-detections, have spatially coincident MIR and 
radio emission, with morphologies that are similar to HRDS Hll regions. Their extrapolated 
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X-band fluxes made them HRDS candidate targets. Their measured X-band continuum 
emission is quite weak, however, so they are probably H II regions that fall below our detection 
threshold. One source, G005. 756— 0.932, i s an IR bubble with coincident NVSS emission. 
It is not included in the IChurchwell et al.l (120071 ) catalog. This source is odd since it has 
the same bubble morphology at 8.0 /im and 24 /xm, whereas all the other Hll regions in 
our sample show 8.0 Atm emission that surrounds, and extends beyond, the 24 /im emission 
(jPeharveng et al.l 120101 ). 



The source G045. 797— 0.347 is probably a region of thermally emitting plasma too 
faint to be detected by the HRDS. It may be a star-forming region that is unrelated to 
the compact sources seen nearby in MIR images. The G017.223+0.397 target is near a 
Lockman et al.l (I19961) "diffuse" Hll region source. It is not, however, a diffuse nebulae: 
Stephenson says th a t it is an evolved star, which is plausible on the basis of our 

non-detection. IStephensoru (119921 ) also claims that G357. 387+0. Oil is a distant, luminous 
star. 

Five other sources seem to be non-thermal emitters. We believe four of them are AGN 
because of the poor alignment of their radio and IR emission. If a source is unresolved at 
~ 1 ' resolution, then the MIR and radio emission should be spatially coincident to within 
~ 30 This is not the case for these sources, which suggests that they are chance alignments 
of AGNs with Galactic MIR sources. The remaining source, G054.096+0.266, is a PWN; it 
is the PSR J1930+1852 pulsar wind nebula associated with SNR 54.1+0.3. Bright at 24 /xm, 
this source is faint at 8.0 /xm, which suggests that line emission may be causing the 24 /im 
brightness. We observed it by mistake. 



4.2. Multiple Velocity Component Sources 

Of the 448 HRDS Hll regions, 129 (nearly 30%) have multiple RRL components in their 
(Hna) spectra at different velocities. We find 105 sources with two components and 23 
with three components. One source has four velocity components. With the 82 " GET beam 
at X-band, it is unlikely that we would detect multiple, physically distinct Hll regions along 
any given line of sight. One of the components must surely come from the HRDS target 
H II region. The other component may stem from more diffuse ionized gas along the line of 
sight. This diffuse component may originate from ionizing photons leaking into the extended 
PDRs produced by nearby, large star-forming regions. 



Even H II regions with well-defined P DRs leak ra diation beyond th e ioniz ation front 
( IZavagno et al.l 120071 : lAnderson et al.ll2010l ). Moreover, lOey fc KennicuttI ( 119971 ) compared 
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the observed Ha luminosity with that predicted from stellar models and concluded that 
up to 50% of the ionizing flux was escaping from Hll regions. The complex morphology 
of a typical HRDS nebula, e.g. Figure IH suggests that within any nebula some zones are 
ionization bounded and others are density bounded. 

If we are to derive the kinematic distances and physical properties of HDRS nebulae, 
we must unambiguously identify their velocities in these multi-component spectra. The 
resolution of this problem is beyond the scope of this paper. Nonetheless, here we make a 
preliminary analysis of HRDS sources located near the W43 star forming H II region complex. 

Over 30% of the HRDS multiple velocity component sources are located within 2° of 
the large, massive star forming region W43. This is a mere ~ 2% of the survey area; it 
represents a significant overabundance of multiple velocity component sources in a very 
small zone of sky. The RRL LSR velocity of W43, measured from a single ~ 3' resolution 
pointing, is 91.6 krns"^ ( L89). This velocity corresponds to a kinematic dis tance of 5.7 kpc 



(lAnderson fc Baniall2009l ). which places W43 at the end of the Galactic bar (iBenjamin et al. 



20051 ) . From studies of extern al Galaxies, we know t hat the end of galactic bars can have 



significant star formation (e.g. iMartin fc Friedlil 119971 ). The entire W43 complex has veloc 



ities r anging from ^ 80 to 105 kms ^, when mapped at high angular resolution with the 



VIA flBalser et al.ll200lh . 



Of the 23 HRDS sources within 0?5 of W43, 21 have multiple velocity components. 
The two single component sources are compact Hll regions that share the same velocity as 
W 43. All the o t her s ources have one velocity component in the RRL velocity range found 



by lBalser et al.l (120011 ). The other component is from the HRDS target source. The origin of 
the multiple- velocity components for the sources in the direction of W43 appears to be low 
density ionized gas along the line of sight. It is at present unclear if this low-density plasma 
is from FUV photons leaking into the extended PDR surrounding W43, from plasma not 
necessarily associated with W43 at the end of the Galactic bar, or from both. 

If the origin of the multiple-velocity components is leaked plasma from W43 itself, we 
would expect the intensity of a source's W43 RRL component to decrease, on average, with 
its distance from the nominal center of W43, {i, b) = (30?78, — 0?03). Here we are assuming 
that the intensity of the RRL emission decreases from the center of W43 outwards, just as 
the continuum emission does. In Figure [5l we plot the intensity of an HRDS source's W43 
RRL component as a function of its distance from the nominal center of W43. Shown is 
the intensity of the RRL component whose velocity is nearest to the 91.6 kms~^ velocity of 
W43. Figure E] shows that these components are on average stronger for sources near W43. 
This is consistent with the hypothesis that these spectral lines are caused by an extended, 
partially ionized PDR surrounding W43. 
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This interpretation, however, is far from secure. Given the distance to W43, an area 1 ° 
in diameter would be ~ 100 pc across. This is extremely large, although not unheard of, for 
a star forming complex. From thermal radio sources such as W43, when RRL emission is 
detected we should also be able to detect radio continuum emission. The lack of extended 
faint continuum emission associated with W43 in the direction of these multiple-velocity 
RRL sources is therefore also a problem for this interpretation. Perhaps we are detecting a 
combination of diffuse emission leaked from W43 and plasma at the end of the Galactic bar. 

It is similarly unclear what is causing the multiple-velocity components for nearly all 
other multiple-velocity HRDS sources. These lines of sight must be individually analyzed in 
order to determine which RRL velocity component is associated with the HRDS compact 
MIR target. This topic will be addressed in future papers. 



4.3. Infrared "Bubble" Sources 



Many Spitzer GLIMPSE "bubble" s ources seem to be associated with Hll regions 
(IBania et al.l l20ld : iDeharveng et al.l |2010| ) . We detected RRLs from the ionized gas in 65 
GLIMPSE "bubble" sources cataloged bv Ichurchwell et al.l J2006I . boOTh . By visually ex- 
amining GLIMPSE images for the entire HRDS zone, 343° < i < 67° with \b \ < 1°, we 
find that 170 of the 420 cataloged GLIMPSE bubbles, ~ 40% of the sample, are associated 
with previously known H II regions. When these are combined with the 65 HRDS detections, 
more than half of all cataloged MIR bubbles turn out to be H II regions with measured RRL 
emission. 

Altogether there are 207 HRDS sources with a bubble morphology; these are targets 
classified as "B," "BB," "PB," or "IB" in Table [H These additional 1 42 HRDS nebulae 
with a bubble morphology that are not in IChurchwell et al.l (l2006l . 120071 ) imply that there 
is a large population of undiscovered bubbles. Within the HRDS zone there are in fact 
hundreds more cataloged and uncataloged GLIMPSE bubbles that have spatially coincident 
radio continuum and MIPSGAL 24 fim emission. They are all Hll region candidate targets. 
They are, however, also very weak continuum sources; all are weaker than the HRDS target 
flux threshold at 9 GHz. 



We speculate that nearly all GLIMPSE bubbles are caused by H II regions (IBania et al. 



20101 ). In the HRDS zone, the combined sample of HRDS and previously known Hll regions 
with measured RRL emission contains ~ 900 nebulae (see §[5]). More than 40% have a 
bubble morphology at 8.0 yum. For longitudes i > 10°, where source confusion is not as 
great as toward the Galactic center, the percentage is > 50%. The bubble morphology is 
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obviously common for Galactic Hll regions. 

It appears that Hll regions with a bubble morphology do not represent any particular 
evolutionary state. Some bubble sources are associated with masers and CS emission, which 
suggests that they are young (see Table [1] and § 14. 5p . For others, their radio continuum 
emission is faint and diffuse, implying an evolved state. The bubble morphology is probably 
a consequence of the density structure of the surrounding ISM. In a turbulent medium, a 
complete bubble would be difficult to form because a density gradient would favor expansion 
in the direction of decreased density. 



On the basis of CO observations, iBeaumont fc Williams! (120101 ) argue that the "bubble" 



sources are in fact physically two-dimensional rings, rather than limb brightened projections 
on the sky of optically thin three dimensional bubble structures. In this scenario, the rings 
are produced by the bipolar flows that occur during star formation. To make a complete 
ring, a bipolar flow source must be oriented at nearly zero inclination angle with respect to 
the line of sight, so the flows of ionized gas are pointing directly toward and away from the 
observer. At other inclination angles these flows make incomplete, bipolar rings on the sky. 

The ring scenario thus predicts two observational consequences that can be tested with 
the HRDS data. First, there should be many more bipolar bubbles than complete bubbles, 
since the latter would be seen only for particular, and rare, inclination angles. Secondly, 
the RRL line widths for bubbles should be larger than those of the Hll region sample at 
large, since both the red- and blue-shifted ionized gas from the bipolar flow contribute to 
the bubble line width. 

The HRDS catalog strongly suggests that bubble sources are three-dimensional objects. 
There are far more nearly complete bubbles, 117 morphological class "B" sources, than 
bipolar nebulae: the HRDS has only 8 morphological class "BB" bipolar bubbles. Bipo- 
lar bubbles are quite rare. Furthermore, the average RRL line width of bubble sources is 
22 kms~^ (see § 15. 3p . which is identical to the average for the entire HRDS sample. Both ring 
scenario predictions are contradicted by the HRDS catalog nebular statistics. We therefore 
conclude that the majority of bubble sources are three-dimensional structures. 



4.4. Notes on Individual Sources 

4.4- i- Possible Planetary Nebulae 

Because PNe are physically small compared to Hll regions, unresolved emission seen 
in GLIMPSE with no nebulosity may be a hallmark of PNe that might be confused with 
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Hll regions. Nearly all of the 8.0 /xm images of HRDS sources show nebulosity. There are, 
however, 13 HRDS RRL sources that appear point-like in the GLIMPSE data. These objects 
may be PNe; they are noted in Table [T] as being PN candidates, "PN?" . These sources are: 
G009. 741+0.842, G012.116+0.076, G012.199-0.034, G016.404-0.575, G017.364+0.519, 
G018.708-0.126, G022.154-0.153, G029.874-0.819, G030.234-0.139, G030.532-0.259, 
G030.663-0.334, G046.017+0.264, and G050.556+0.04 5. Four of these sources were 
detected as part of the RMS survey (lUrquhart et al.l 120091 ) and these authors classify them 
as PNe. Since all targets classified as point sources ( "PS" in Table [T]) are similar in this 
regard, it seems likely that most of them are also PNe. 

Based on their 8.0 /im morphology, four additional sources are also likely to be PNe: 
G003.449-0.647, G016.228-0.368, G017.414+0.377, and G026.317-0.01. Two of these 
sources, G017.414+0.377 and G026. 317— 0.012, are observed as faint rings at 8.0 fim. Their 
low 8.0 fim intensity and morphology separate them from the other HRDS 8.0 fim rings. For 
G026. 317— 0.012, a central point source is detected within the ring. The two other sources, 
G003.449-0.647 and G016.228-0.368, a re not detecte d at 8 .0 /xm. They have a "ring" 
morphology at 24 /im and are listed in the iMizuno et al.l (120101 ) catalog of MIPSGAL ring 
sources. These authors suggest that a large percentage of sources observed as rings at 24 fim 
are PNe or evolved stars. 



The HRDS source G031. 727+0.698 is also hsted in the iMizuno et al.l tom catalog. 
This source has a negative LSR velocity of —39.2 kms~^ and it is therefore located in the 
outer Galactic disk. Its large distance gives it a high intrinsic luminosity and a large physical 
size, which means that it cannot be a PN. Th e minimal over l ap be tween the HRDS, which 
contains 207 8.0/im rings (bubbles), and the IMizuno et al.l ( 120101 ) catalog of 24 /im rings 
implies that while 8.0 /im rings nearly always surround Hll regions, the rings detected at 
24 /im are associated with different objects. 

Due to their expansion and higher el ectron temperature, the line widths o f PNe are gen- 
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(Il997l ). for example, find an average hydrogen recombination line width of 44.1+5.8 km s~^ for 
a sample of six PNe. The 17 HRDS sources identified here as being possible PNe have an av- 
erage line width of 25 kms~^, which is near the mean of the entire HRDS sample, 22 kms~^. 
Three sources, however, do have broad RRLs: G026.317-0.012, G029.874-0.819, and 
G050. 556+0. 045 have line widths of 38.0, 40.6, and 48.8 kms~^, respectively. Surpris- 
ingly, six HRDS sources listed as possible PNe have narrow line widths, < 17kms~^: 
G009. 741+0.842, G016.404-0.575, G017.364+0.519, G017.414+0.377, G018. 708-0.126, 
and G022. 154— 0.153. More observations are needed to better understand the nature of 
these sources. 
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4.4-2. Non-detections at 8.0 fim 

Eight sources have 24 /iin emission but no detectable 8.0 fim emission: G003. 449— 0.647, 
GOll. 573+0.340, G016.228-0.368, G019. 786+0.285, G023.513-0.244, G023.849-0.196, 
G029.019+0.170, and G049.507-0.520. These objects are hsted in Table [U as a non- 
detection, "ND". That we detect RRL emission from these targets points to a thermal 
origin. Two of these sources, G003. 449— 0.647 and G016. 228— 0.368, are identified above as 
being PNe candidates. Finally, G049.507-0.520 has no associated GLIMPSE or MIPSGAL 
emission but is likely associated with W51 (see below). The four remaining sources all have 
an extended, curved morphology at 24 fim and are thus unlikely to be PNe. 



4.4- 3. Nuclear Disk 



We found nine sources whose velocity < —200 kms ^ and Galactic location of £ ~ 
359 ° place them in the nuclear disk of the Milky Way: G358.530+0.056, G358.552-0.025, 
G358.616-0.076, G358.652-0.078, G358.680-0.087, G358.694-0.075, G358.720+0.011, G358.827+0.085, 
and G359. 159-0.038. Before the HRDS, only three Hll regions were known to share this loca- 
tion and velocity range: G358.623-0.066 at -212 kms'^ G358. 797+0.058 at -20 6.6 kms-\ 



and G358. 974— 0.021, which has two RRL components at —193.3 and —5.4 km s ^ ( Caswell fc Haynes 
19871 : iLockman et al.l 119961 ). The location observed by ICaswell fc HaynesI (Il987l ). {i, b) = 
(358.623, -0.066), lies in between the two HRDS sources G358.652-0.078 and G358.616-0.076; 
it is not a distinct Hll region. These HRDS sources are distributed within a Galactic zone 
~ 30' in diameter, or ~ 75 pc assuming a distance of 8.5 kpc. Most of these nebulae have a 
similar compact morphology (we classify five of them as "Compact" in Table [1]) suggestive 
of relatively early evolutionary stages. 



4.4.4. G036. 914+0.489 

We re-observed G036. 914+0. 489 in order to resolve confusion about its velocity. L89 
found two RRL velocities for his G036.914+0.489 target: -30.3 and +47.0 kms"\ There 
are in fact two distinct radio and IR peaks associated with this region: G036. 870+0. 462 
and G036. 914+0. 485. The velocities of both HRDS sources are near —30 kms~^, implying 
a large distance from the Sun of ~ 17 kpc (Anderson et al., 2011b, in preparation). We do 
not detect a second emission component from either source. 
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4-4-5- Coincidence with Infrared Dark Clouds 



Some HRDS nebulae coincide with the locations of IR dark clouds (IRDCs). IRDCs are 



assoc iated with the earliest phases of star formation ( IChambers et al.ll2009t iRathborne et al. 
20101 ) ■ If the nebulae and IRDCs are truly spatially co-incident, then these regions are 



candidates for the youngest sources in the HRDS. The best examples of this phenomenon 
are G009.875-0.749, G016.361-0.209, G018.832-0.300, G027.366-0.164, G028.394+0.076, 
G029. 770+0.219, G030.345+0.092, and G030.378+0.106, although this is not a complete 
list. These HRDS sources generally have a small angular size and are frequently classified 
as "compact." This suggests a small physical size and an early evolutionary state. The 
physical association between these HRDS nebulae and IRDCs, however, is not always clear. 
For example, G028. 394+0. 076 is entirely surrounded by a very large IRDC. We detect two 
velocities for this source: 42.8 an d 86.2 kms~^. The 86.2 kms~^ velocity is similar to that of 
the IRDC, which is 75.8 kms~^ (jsimon et al. 2006). The more than 10 kms~^ difference is 
surprising, however, if the two are indeed associated. The nature of the relationship between 
this H II region and the IRDC remains unclear. 



4-4-6- Strange Line Profiles 



Two sources have double velocity component line profiles that are different from all other 
HRDS targets. G028.764+0.281 has hydrogen RRLs at 104.5 and 87.2 km s'^ G038. 120-0.227 
has its RRLs at 54.7 and 89.3 km s~^. Both sources have line profiles with a narrow com- 
ponent of line width < 13 kms~^ and a broader component. For G028. 764+0. 281, the 
narrow component is at 104.5 kms~^ and has a FWHM line width of 12.5 kms~^. For 
G038. 120—0.227, the n arrow componen t has a velocity of 54.7 kms~^ with a FWHM line 



width of 12.6 kms ^. iBronfman et al.l (119961 ) measured a CS velocity of 83.5 kms ^ for 



G038. 120—0.227, which implies that the narrow 54.7 km s~^ line is from another source 
along the line of sight. In both cases it is not obvious w hat additional sour ce along the line 
of sight is producing such a narrow recombination line. lAdler et al.l (Il996[ ) argue, however, 
that the narrow RRL lines that they see in W 3 are caused by partially ionized hydrogen. 



^.^.7. Unusual Morphology 

One source, G349. 579— 0.680, has a morphology that is not seen in any other HRDS 
target. Its compact GLIMPSE 8.0 /im emission has the appearance of a bright-rimmed 
cloud, whereas its MIPSGAL 24 /xm emission is more extended and diffuse. G349. 579— 0.680 
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is near the center of a large (nearly 1 ° in diameter) 8.0 fim bubble located at {i, b) = (349 ?5, 
— 0?6). There is no published RRL velocity for the large bubble although it appears to be 
an Hll region based on the mid-IR morphology. The —19.4 kms~^ velocity we measure for 
G349.579-0.680 may be that of the large bubble. 



Confusion with Other Hii Regions 

In three cases an HRDS target may be part of a larger, more evolved Hll region. For 
example, G347.536+0.084 is part of the ionization front of the large Hll region centered 



at (£, b) = (347?5, +0?2). ICaswell fc HaynesI (Il987l ) found RRL emission in two locations 



along the ionization front at velocities of —97 and —96 kms ^. We measure a velocity of 
-102.5 for G347.536+0.084. One region, G049.507-0.520, has no GLIMPSE or MIPSGAL 
counterpart; it is a region of radio continuum emission near W51 {i, b) = (49?5,— 0?4). 
Its 63.6 kms~^ RRL velocity is near that of W51, for which L89 measured a velocity of 
58.2 kms~^ Finally, G006.014-0.364 may be a second ionization front for G005.899-0.427, 
for which L89 measured a velocity of 10.1 kms~^. We measure a velocity for G006. 014— 0.364 
of 14.2 kms-^ 



4.5. Astrophysical Context 

Many of the HRDS sources appear in catalogs made from existing Galactic surveys, but 
not as H II regions. The 448 nebulae detected here in RRL emission for the first time, when 
correlated with previous catalogs made at multiple wavelengths, can provide insight into the 
properties of star forming regions. Here we focus on the astrophysical implications of HRDS 
associations with methanol masers and dense molecular clouds traced by CS emission. 



4-5.1. Methanol Masers 



Hll regions are expected to be associate d with methanol masers in the earliest phases 
of high mass star formation (lEllingsenl l2006l ) . Although many of our HRDS nebulae have 
been observed in meth anol maser transitions, there are relatively few detections. Using the 
Pestalozzi et al.l (120051 ) compilation of 14 methanol maser studies, we find that only 10% 
(46/448) of the HRDS sources have detected methanol maser emission within 2' of their 
position. Table [1] lists the velocity for the strongest maser component. This low correlation 
with methanol masers suggests that either our sample of Hll regions is on average not 
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extremely young o r that the maser emission is too weak to be detectable. For comparison, 



Walsh et al.l ( Il998l ) found that 38% (201/535) of the sources in their UC Hll region sample 



have methanol maser emission. This suggests that the HRDS nebulae are not as young as 
their UC sample. We defer a more detailed study of the evolutionary state of the HRDS 
nebulae to a future publication. 



4-5.2. Dense Gas 



Tracers of d ense molecular gas, su ch as CS, are hallmark chara cteristics of young 



Hll regions (e.g., IChurchwell et al.l Il990f ) . The [Bronfman et al.l (Il996[ ) CS (2 — )■ 1) sam- 



ple contains 110 of o ur HRDS targets. Bronfman et al. ( 19961 ) observed all IRAS point 
sources satisfying the Wood fc ChurchwenT i 1989a ) IRAS color criteria for UC Hll regions. 
Roughly half (57/110) of these HRDS sources were detected in CS. Overall, the veloc- 
ities are in good agreement; the mean difference between the CS and RRL velocity is 
4.2 ± 5.9 kms~^. There are, however, four sources with velocities that differ by more than 
10 kms-i; G023.585+0.029, G028.304-0.390, G349.437+1.058, and GOOO.382+0.017 . 
Furthermore, there are two double-velocity RRL sources whose CS velocity lies in between 
the RRL velocities: G352.521-0.144 and G049.998-0.125 . Excluding these sources, the 
average velocity difference is 3.0 ± 2.3kms~^. Thus, with some exceptions, the CS and RRL 
velocities are in very good agreement. 

The CS non-detections are equally interesting since they imply that the Hll region 
is evolved and has displaced or dissociat ed its dense mol e cular gas. IRAS point sources 
satisfying the color criteria employed by iBronfman et al.l (119961 ) have generally been as- 
sumed to be UC H l l regio ns; in many cases this assumption has been proven valid (e.g. 
Wood &: Churchwelll Il989bl ). That half of the HRDS sources targeted were not detected 
in CS, however, implies instead that many of the HRDS nebulae are distant and evolved. 
The CS sensitivity limit may have precluded the detection of faint CS emission from distant 
nebulae. 



5. Properties of the HRDS Nebulae 

We derive here the statistical properties of the HRDS sources. Although the HRDS 
found 448 Hll regions, the number of physically distinct nebulae represented by the 603 
discrete hydrogen RRL components is not well known. Thermal radio sources can often be 
resolved into several apparently physically distinct emission regions, each having a somewhat 
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different position and RRL velocity. Furthermore, the most massive star-forming complexes, 
W43 for example, are extended and fragmented into many sub-clumps of localized star forma- 
tion, which together can ionize a very large zone. Many of our multiple velocity component 
targets may be detecting such low density ionized gas in addition to RRL emission from 
another, physically distinct nebula (see § 14. 2p . Nonetheless, here we follow the convention 
established by L89 and assume that each of the HRDS 603 RRLs is produced by a distinct 
object. 

We also wish to compare the properties of the newly discovered HDRS nebulae with 
those of the previously known Hll region census. There is, however, no extant Hll region 



catalog that is suitable for this purpose. We ha ve therefore compi 



known H II regions frorn existing RRL surveys ( Reifenstein et al 



Caswell fc Havne^ 



20031: 



Sewilo et al. 



987: Lockman 1989: Lockman et al 



1996 



ed a catalog of previously 



1970: Downes et al 



Araya et al 



2002 



1980 



Watson et al. 



20041 ). This catalog improves upon previous efforts not only because it 



uses the most recent observations but also because it accounts for duplicate sources. 

Many Hll regions appear in multiple catalogs, so extreme care must be taken when 
combining them into a single compilation. We visually inspect the radio continuum and 
MIR emission from all previous RRL observations of HIT regions to help identify sources that 
were observed by multiple authors. If observed positions are part of the same contiguous 
radio continuum and MIR zone of emission and have RRL velocities within 5 kms~^ of one 
another, we deem them to be observations of the same object. For such matches we only 
keep the source data from the most recent catalog, because the newest recombination line 
data are almost invariably the most accurate. We have also removed known SNRs, LBVs, 
and PNe from the sample through correlation with the SIMBAD database using the nominal 
HRDS target positions and a 5 ' search radius. 

We refer hereafter to this catalog of H II regions previously known to reside within the 
HRDS longitude and latitude range as the "Known" sample. For each source, the Known 
sample catalog contains the Galactic longitude and latitude, the RRL LSR velocity, and 
the FWHM line width. The Known sample data can be found online as an ASCII text file 
1^; Anderson et al. (2011a, in preparation) give the full details. There are 456 unique Hll 
regions in the Known sample, although there is likely some residual contamination. With its 
448 unique positions, the HRDS has thus doubled the sample of H II regions in the survey 
zone. 



^http://go.nrao.cdu/hrds. This catalog is as yet a work in progress. It still is likely to contain some 
residual contamination and duplicate entries. 
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5.1. Galactic Distribution 

The Galactic longitude distribution of the HRDS and Known samples, shown in Figure|6l 
are very similar. In Figure |6l the HRDS sample is shown as dark gray and the Known sample 
is shown as a dotted outline. A Kolmogorov-Smirnov (K-S) test shows that the distributions 
are not statistically distinct. The HRDS Galactic latitude distribution, shown in Figure [TJ 
peaks closer to the Galactic mid-plane compared with the Known sample. A Gaussian fit to 
the two distributions shows that both the HRDS and the Known sample have essentially the 
same FWHM scale-height of 0?6. A K-S test of the two samples, however, reveals significant 
statistical differences, entirely due to the dif ference in peak la titude. It is not clear why the 



two distributions peak at different latitudes. lAndersonI (120101 ) found that the nebulae in our 



pilot study for the HRDS were on average more distant than the Known sample Hll regions. 
The HRDS may simply be sampling a different area of the Galaxy. 



5.2. Galactic Structure 

The fundamental result of any RRL survey is the LSR velocity of a nebula. Thus the 
fundamental map of the observed data is its (£, v) diagram. The HRDS and the previously 
known H II regions together give a census of the Inner Galaxy, the zone within 70 ° of the 
direction to the Galactic center, that shows clear evidence for an ordered pattern of Galactic 



structure in {i,v) space (iBania et al.ll2010l ). This was asserted by previous researchers, but 
our new census, which has a factor of two more objects, now unambiguously shows the kine- 
matic signatures of the spatial distribution of H II regions, a concentration of nebulae at the 
end of the Galactic bar, and nebulae located on the kinematic locus of the 3 kpc Arm. The 
HII region Galactocentric radial distribution shows very narrow (A Rg ^ 1 kpc) concentra- 
tions of star formation at Galactic radii of 4.25 and 6.00 kpc. This was k nown before but th e 



HRDS sharpens the contrast and now makes this structure very robust (IBania et al.ll2010l ). 
More detailed Galactic structure studies that use kinematic distances for HRDS nebulae 
determined by Anderson et al. (2011b, in preparation) are underway (Bania et al. 2011, in 
preparation) . 

The HRDS has 34 negative- velocity H II regions in the first Galactic quadrant at i > 10 °. 
This represents 7% of all HRDS nebulae in this zone. First-quadrant sources with negative 
velocities are unaffected by the kinematic distance ambiguity; assuming circular rotation 
about the Galactic center, a negative RRL velocity implies that a nebula is beyond the Solar 
orbit at a Galactocentric radius of more than 8.5 kpc. Before this survey, there were only 
seven known H II regions with negative- velocities in this region. We have therefore increased 
the sample of negative- velocity Hll regions thought to lie beyond the Solar orbit in this 



-27- 



Galactic zone by almost 500%. The negative-velocity sources are fainter on average when 
compared with the HRDS sample as a whole. Their average continuum intensity is ~ 250 mK, 
whereas the HRDS average is ~ 300 mK. One negative-velocity source, G032. 928+0. 607, has 
a continuum intensity of ~ 700 mK, making it one of the brightest HRDS nebulae. Because 
of their extreme distances, these sources are quite luminous. They lie along the (£, v) locus of 
the Milky Way's "Outer Arm" (IBania et al.ll2010l ). These regions of massive star formation 
in the outer Galaxy can provide important constraints on Galactic chemical evolution. 

The HRDS has nine sources in the Galactic center direction with RRL velocities smaller 
than —150 kms~^. This suggests that they lie in the "Nuclear Disk," the region of high 
velocity gas in the inner ~ 0.5 kpc of the Milky W ay. These sources all lie between 358 ?5 < 
i < 359 ?2 and \b \ < 0?1. iLockman et al.l (119961 ) found two "diffuse" nebulae in this part 
of the sky, which are the only previously known nebulae with similar velocities. We find 
only one source, GOOO. 729— 0.103, with a velocity, 105.3 kms~^, that could place it in the 
red-shifted side of the nucle ar disk. Although th ese ten sources show that star formation is 
ongoing in the nuclear disk, ICaswell et al.l ( l2010l ) find no methanol masers with nuclear disk 
velocities. 



5.3. Radio Recombination Line Widths 

The ~ 25kms~^ line widths of Hll region RRLs are set by a combination of thermal, 
turbulent, and ordered motions (e.g., flows). This line width greatly exceeds the thermal 
width expected from the ~ 10^ K temperature that is typical of H II region plasmas. H II 
region RRLs are thus significantly broadened by turbulent and ordered motions. In Figure |H1 
we show the RRL full width at half maximum (FWHM) line width distribution for the HRDS 
sample (dark gray) and the subset of the L89 sample from the HRDS Galactic longitude and 
latitude range (dotted outline). We compare the HRDS with a single sample, that of L89, 
rather than the Known sample in order to minimize any possible systematic effects stemming 
from different observing and data analysis techniques. 

The typical FWHM line widths for the HRDS nebulae are narrower than those of the 
L89 Hll regions. The distribution for the HRDS sample peaks at ~ 20 kms~^, whereas the 
L89 sample is noticeably broader and shows peaks at 23 and 26 kms^^. The formal means 
and standard deviations are 22.3 ± 5.3 for the HRDS compared with 26.4 ±8.1 for L89. A 
KS test of the HRDS line widths (603 sources) and the L89 line widths (462 sources) shows 
that the two distributions are statistically distinct to a high degree of certainty. Some of 
the difference between the two distributions may be due to the different spectral resolutions, 
which are 1.86 and 4 kms~^, respectively, for the HRDS and L89 surveys. The L89 line 
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widths could also be broader because the 140 Foot telescope's bigger beam (195" compared 
with the GBT's 82 " ) sampled a much larger volume of gas. 

The HRDS has 35 nebulae with extremely narrow, < 15 kms~^, line widths. These 
line widths can be used to estimate the maximum electron temperature of the Hll region 
plasma. Since turbulence and ordered motions only broaden the RRLs, interpreting the line 
width as being purely thermal sets a robust upper limit on the nebular electron temperature. 
Extremely narrow line widths, < 15kms~^, ir nply "cool" nebulae with electron temperatures 
< 5,000K. Such nebulae were discovered by IShaver et al.l (119791 ). They were first detected 



in appreciable numbers by L89, who speculated that cool nebulae may be common and that 
few were known because of the sensitivity limits of the existing RRL surveys. The HRDS, 
however, shows that cool nebulae are indeed rare in the Galaxy: only 6% of HRDS sources 
have line widths smaller than 15 kms~^. This fraction is, however, three times higher than 
that found for the L89 sample. 

We detect only three nebulae with large, > 40 kms~^, line widths. Ignoring any contri- 
bution from ordered motions and assuming that the thermal and turbulent components add 
in quadrature, lines broader than 40 kms~^ imply turbulent motions of order 35 kms~^ for 
a typical 10, 000 K Hll region plasma. Some broad line width sources, however, may in fact 
be PNe (see UM- 

As the spectral signal-to-noise decreases, the line width derived from Gaussian fitting 
becomes increasingly uncertain. Furthermore, with low signal-to-noise, double-velocity com- 
ponents may be not be well-separated and an erroneously large line width might be derived. 
In Figure [H we plot the line width as a function of the line intensity for the HRDS (left 
panel) and the L89 samples (right panel). Figure [9] shows that the line width distribution 
for low intensity HRDS sources is similar to that of the entire distribution. This is not the 
case for the L89 sample — many broad line sources are of low intensity. We conclude that at 
centimeter- wavelengths Galactic Hll regions have line widths < 35 kms~^. The HRDS has 
only eight lines (1%) broader than 35 kms~^, and only four lines broader than 36 kms~^. 
Line widths greater than 35 kms~^ should be regarded with suspicion because they may be 
the result of a low signal-to-noise detection, represent blended velocity components, or a sign 
that the source is a PN and not an H II region. 



5.4. Continuum Angular Size 



The distribution of the continuum angular sizes for HRDS nebulae is shown in Figure [TO] 
which plots the geometric mean of Gaussian fits to the RA and Dec continuum scans. Al- 
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together, we have a useable continuum measurement in either RA or Dec. for 441 HRDS 
nebulae. Here we use the fitted FWHM angular sizes uncorrected for the GET beam. Typi- 
cal errors for the fits are 5 to 10 ". Assuming Gaussian source and beam shapes, the observed 
and true angular sizes are related by the expression: 6"^ = 6"^ + Of , where 6 is the observed 
size (shown in Figure [TOj) . 9s is the true source size, and 9b is the beam size. The vertical 
dashed line in Figure [TO] flags the 116 " size that would be measured by the GET for a source 
that just fills its X-band beam. After correcting for the beam size, the majority of the HRDS 
nebulae, 235 of 441 (53%), are unresolved by the 82" GET beam. Eleven sources (~ 5%) 
have angular sizes larger than 250 "; these are not shown in Figure [TOl 

Although small in angular size, many HRDS nebulae lie at extreme distances from the 



Sun (jAndersonll2010l . Anderson et al., 2011b, in preparation). Thus many of these angularly 
small nebulae are actually large in physical size. If the HRDS nebulae have on average large 
physical sizes, this would imply an evolved state and would explain why the HRDS targets 
are in general not correlated with methanol masers and CS emission (see 14. 5p . 



5.5. Planetary Nebula Contamination 

There is likely some contamination from PNe in the HRDS sample because our selection 
criteria locate thermally emitting sources and not just Hll regions. We have, however, 
identified only 17 HRDS sources that might be PNe (see § 14.4. ip . Thus fe wer than 4% 



(17/44 8) of the HRDS sources are potential PNe candidates, which confirms lEania et al. 



(I2OIOI ) who argued that the level of contamination in the HRDS by PNe must be quite 
small. This conclusion was based on considerations of Galactic structure, scale height, RRL 
line widths, and the RRL line-to-continuum ratios (i.e. nebulae electron temperatures). In 
particular, the clear evidence for Galactic structure seen in the [i, v) and Galactic radial 
distributions of HRDS nebulae implies that the level of PNe contamination in the HRDS 
sample must be minimal. Eecause PNe are an old stellar population their Galactic orbits are 
well-mixed. PNe show, therefore, no structure in their Galactocentric radial distribution and 
their Galactic (£, v) distribution is a scatter plot constrained only by velocities permitted by 
Galactic rotation. Any PNe contamination of the HRDS sample must therefore be very small, 
otherwise these interlopers would suppress the unambiguous signal of Galactic structure seen 
in the HRDS Galactocentric radial and {£, v) distributions. 
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5.6. Survey Completeness 



By assuming a flux distribution, one can assess the completeness of a flux-limited survey. 
There are good reasons to expect a power law for the HRDS fluxes. For example, power 
law flux distributions have been found for the MSX, GLIMPSE, and IRAS point source 
catalogs. A power law distribution of fluxes reflects an underlying power law for source 
luminosities together with a relatively structureless pattern of source locations across the 
Galaxy. For H II regions, it is well-known that a power law is a good approx imation to the 
luminosity function (e.g., ISmith fc Kennicuttlll989l : iMcKee fc Williamslll997l ) and that they 
are relatively smoothly located across the Gal axy. While there are strong peaks in the H l l 
region Galactocentric distance distribution (cf. lAnderson fc Baniall2009l : iBania et al.l 120101 ). 
their heliocentric pattern shows no structure. 

We attempt to assess the completeness of the HRDS by evaluating the integrated con- 
tinuum flux distribution calculated using Equation [T] and assuming a power-law distribution 
of fluxes. 

Because the radio continuum data suffer from confusion and possible zero point uncer- 
tainties, we also use the RRL data to estimate the continuum flux of the HRDS nebulae. 
To do this we assume an X-band l ine-to-continuum ratio of 0.09. This ratio is the average 
value of the iQuireza et al.l (l2006al ) sample of Galactic H II regions, derived using their 28 
highest quality measurements (their quality factor values of "A" or "B" for both the line 
and the continuum data). We then use Equation [1] and the fitted RA and Dec continuum 
scan widths to get the integrated continuum fluxes. In both cases we are assuming that the 
sources are Gaussian in both RA and Dec and well-characterized by the cross scans. We 
are also implicitly assuming that the HRDS nebulae do not suffer greatly from optical depth 
effects (see gSJ^D- 

Histograms of the integrated source flux distribution for HRDS nebulae are plotted in 
Figure [TTl where the bottom panel shows the fluxes derived using the continuum data and 
the top panel shows continuum fluxes estimated from the RRL data. Solid lines are power- 
law fits to a panel's data for the range 200 to 1,000 mJy The dashed line shows the fit to 
the other panel's data. We estimate the completeness limit using the histogram bin where 
the fit deviates significantly from the fitted power law; these limits are shown as vertical 
dotted lines. The fit to the continuum flux distribution has a slope of —0.8 and shows a 
completeness limit of ~ 280 mJy This fit is, however, rather poor (reduced x^=15). The fit 
quality of the RRL estimated continuum flux is much better (reduced x^=3). It is best-fit 
with a slope of —1.0 and shows a completene ss limit of ~ 180 mJy The HRDS completeness 
limit of 120 mJy cited by iBania et al.l (120101 ) was based on an early analysis of sources from 
i = 30° to 50°. It appears that the completeness limit is lower in this region of the Galaxy, 
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possibly because these candidate targets were identified using the more-sensitive VGPS data. 

There is a large difference between the two distributions shown in Figure [TTJ Some of 
the discrepancy is caused by assuming a single line-to-continuum ratio for the entire sample. 
In reality the line-to-continuum ratio varies between Hll regions. Perhaps a larger issue, 
however, is that many of our targets are in locations with significant continuum emission 
structure. It is often difficult to separate the target's emission from that of other sources 
and/or noise due to sky fluctuations. In many cases it is also difficult to establish a reliable 
zero point for the continuum intensity. Because of these concerns, even though we must 
assume a line-to-continuum ratio, using the line intensity as a proxy may in fact be a better 
estimate of the actual continuum flux of faint sources. 

The completeness limit derived here is over twice that expected from our HRDS target 
list. This discrepancy almost certainly stems from the uncertainty in the extrapolation of 
20 cm radio continuum fluxes to predicted 3 cm (X-band) GBT fluxes. We found empirically 
that the flux density actually measured at X-band can differ by a factor of two from the 
extrapolated flux density. 



5.6.1. Spectral Type Completeness 



At X-Band, RRL emission propagates on average through an optically thin medium 
and thus we can detect Hll regions on the far side of the Galaxy. Here we estimate the 
distance to which we can detect H II regions ionized by stars of various spectral types given 
the derived HRDS completeness limit of 180 mJy. To do this we first must calculate the 
luminosity dependence of H II regions for a range of ionizing spectral types. 

The spectral type of the exciting star of an H II re gion, assuming ionization by a single 
star, can be estimated from its radio flux (IRubinl 119681 ): 



iViy ^ 4.76 X 10"» I ^ 1 ( 



48 



Jy 



K 



-0.45 



(2) 



where A'ly is the ionization rate, the number of Lyman Continuum ionizing photons emitted 
per second. Si, is the nebular radio flux density, and Tg is the electron temperature. This 
equation is almost unaffected by source geometry and variations in electron density within 
the source, but does assume that the continuum emission is optically thin. 

Using the Lyr nan continuum photon emission rates for stars of spectral type 03 through 
BO. 5 calculated by lSternberg et al.l (120031 ) . we convert the Equation [2] A^iy values into ionizing 
star spectral types. We plot in Figure [12] the expected flux density of H II regions ionized 
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by single stars of spectral type 03 to BO as a function of the distance to the Hll region. 
The gray curves in Figure [12] show the estimated flux for the range of electron temperatures 
found in H II regions, ~ 5, 000 to ~ 10, 000 K. The dashed line shows the 70 mJy extrapolated 
9 GHz flux limit of the HRDS target list. The dash-dot line shows the actual HRDS 180 mJy 
flux limit derived in §5.61 The HRDS is complete for all Hll regions surrounding single 
0-stars within 12kpc of the Sun (dotted line). Being flux-limited, however, the HRDS is 
complete to different distances for each spectral type. The HRDS can detect Hll regions 
ionized by single 07 stars across the entire Galactic disk. Since the ionizing flux decreases 
rapidly for later spectral types, however, we can only detect Hll regions ionized by Bl stars 
at a distance of only a few kpc. 



5.6.2. Optical Depth Effects and Confusion 

Young Hll regions may be optically thick at 20 cm, an effect which could affect our 
completeness. We use the 20 cm VGPS and NVSS data to estimate the X-band flux. If a 
nebula is optically thick at 20 cm then its extrapolated X-band flux would be underestimated. 
The source might then fall below our flux threshold and so not make it onto the HRDS target 
list. The flux density of an Hll region rises as z/^ until it peaks at a frequency, z/q, called the 
"turnover frequency." At frequencies less than z/q the nebula is optically thick whereas at 
frequencies greater than it is optically thin. In the optically thin regime, the flux density 
has a very weak dependence on frequency, S oc z/~°'^. The turnover frequency is: 

T,\-»"V EM ^»"™ 



= 0.3045 ^^j IGHzl. (3) 

where EM is the emission measure and is the electron temperature. For compact Hll 
regions with Tg ~ 10^ K and EM 10^ cm~^ pc, the turnover frequency is near 1 GHz. Some 



UC a nd hyper-compact H II regions may have emission measures in excess of 10^ (iKeto et al. 



20081 ). which leads to turnover frequencies near 15 GHz. Thus, optical depth effects become 
important for young UC H II regions with high emission measures and for observations at low 
frequencies. Our survey therefore is less complete for the youngest UC and hyper-compact 
Hll regions that are optically thick at 1.4 GHz because we have not targeted these nebulae 
in our observations. 

While we make no attempt to prove empirically that our target location method flnds 
all possible H II region candidates, it does not seem to be severely affected by non-detections 
in the radio and IR surveys or by source confusion. In the HRDS survey zone we flnd over 
1,000 Hll region candidate targets that have spatially coincident MIR and radio continuum 
emission with a similar morphology and angular extent. About 40% of these objects have 
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20 cm flux densities below the 85 mJy HRDS 20 cm flux density cutoff. Our method thus 
easily finds a large number of sources fainter than the fiux threshold wc used for the HRDS. 
We conclude that the apparent brightness of sources in the radio and IR data does not 
severely impact our completeness. 

Confusion may also limit the completeness of this survey. Because of the large number 
of sources found near W43 at {I, b) = (30?78,-0?03), which is the most complicated emission 
region within the survey zone, we do not believe confusion significantly affects our ability to 
visually identify Hll region candidates. In large part this is because that radio emission is 
relatively unconfused compared to the IR emission. The most significant problem affecting 
our completeness is likely the errors in the X-band ffuxes extrapolated from the 20 cm radio 
continuum data because this has a strong impact on which sources were observed. 

6. Summary 

The Green Bank Telescope H II region discovery survey at 9 GHz has doubled the number 
of known Hll regions in the Galactic zone 343° < i < 67° with | 5 | < 1°. The HRDS targets 
were selected based on spatially coincident 24 /xm and 20 cm emission of a similar angular 
extent. Our RRL detection rate of 95% proves that this criterion is extremely effective 
in identifying Hll regions. The derived survey completeness limit, 180 mJy at X-band, is 
sufficient to detect all optically thin Hll regions ionized by single 0-stars at a distance of 
12 kpc from the Sun. These recently discovered nebulae share the same distribution on the 
sky as does the previously known census of Galactic H II regions. On average, however, the 
new nebulae have fainter continuum ffuxes, smaller continuum angular sizes, fainter RRL 
intensities, and smaller RRL line widths. 

Nearly 30% of the HRDS Hll regions (129 nebulae) have multiple hydrogen RRLs in 
their spectra at different velocities. Moreover, over 30% of these multiple velocity component 
sources are located with 2°of the large, massive star forming region W43. Low density ionized 
gas produced by FUV photons leaking into the extended PDR surrounding W43 appears to 
be responsible for nearly all the multiple velocity components detected in this region of 
sky. We speculate that the majority of the HRDS multiple velocity component RRL spectra 
stem from similar low density PDR plasmas distributed over a large region of star formation. 
These zones are so large in angle that they overlap the fines of sight to our HRDS targets. 

The HRDS has 207 Hll regions with a bubble morphology at 8.0 /xm in the GLIMPSE 
survey, 65 of which have been cataloged previously. By comparing the positions of bubbles 
detected in GLIMPSE with the Known sample of H II regions, we find that about half of all 



- 34 - 



Galactic Hll regions have a bubble morpholog y in the infrared. We speculate that nearly all 



GLIMPSE bubbles are caused by Hll regions ( iBania et al.ll2010l ). The bubble morphology is 



probably not associated with any particular evolutionary stage, but is rather an indication 
of the homogeneity of the ambient environments of these H II regions. The HRDS statistics 
on the distribution of different morphologies and the line widths of bubble sources strongly 
suggest that these objects are limb brightened projections on the sky of optically thin three 
dimensional structures, rather than two-dimensional rings produced by bipolar flows. 

The HRDS has 34 Hll regions with negative LSR velocities. These first Galactic quad- 
rant nebulae he beyond the Solar orbit at (Iq > 12 kpc in the outer Galactic disk, Rq ^ 
9kpc, placing them in the high-z, warped Outer Arm. Previously there was only a single 
Hll region known in this zone. This new sample of outer Galactic disk Hll regions paves 
the way for future studies. Because metals are the main coolants in the photo-ionized gas, 
the H II region heavy element abundance can be derived by proxy using the nebular electron 
temperature, Tg. The electron temperature of Hll regions is directly related to the amount of 
processing of the gas and dust from which it formed, and thus to Galactic chemical evolution 
(GCE). There are relatively few Hll regions with accurately derived valu es, especially fo r 



the Rg ^ 10 kpc region which is critical for constraining models of GCE ( iFu et al.ll2009l ). 
Determining electron temperatures for these Outer Arm and other HRDS nebulae will be of 
great importance for future GCE studies. 

The HRDS also has ten Hll regions whose position and velocity place them in the 
nuclear disk. Previously there were only two Hll regions known to lie in this part of the 
Galaxy. These HRDS nebulae show conclusively that active star formation is occurring in 
the innermost part of the Milky Way. 

The distribution of HRDS nebular line widths peaks at ~ 20kms~^, which is smaller 
than the ~ 25kms~^ line widths found for the previously known Hll region sample. The 
combined census distribution shows that Galactic H II regions at 9 GHz have line widths 
smaller than 35kms~^. Line widths greater than this may be the result of multiple com- 
ponents along the line of sight, a low signal-to-noise detection, or misidentification of a 
planetary nebula as an Hll region. The HRDS has 35 nebulae with extremely narrow, 
< 15kms~^, line widths. These narrow line widths imply "cool" nebulae with elec tron tem- 



peratu res < 5, 000 K. Such nebulae were first detected in appreciable numbers by iLockman 



( 119891 ). who speculated that cool nebulae may be common and few were known because of 
survey sensitivity limits. The HRDS shows, however, that cool nebulae are indeed rare in 
the Galaxy: only 6% of HRDS Hll regions have line widths smaller than 15kms~^. 

Finally, the majority of the HRDS Hll regions are unresolved with the 82" GET beam 
at X-band. Despite their small angular size, however, it would be wrong to assume that these 
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nebulae are in an early stage of their evolution. The HRDS sample shows a poor correlation 
with dense gas and methanol masers, which are usually associated with young H II regions. 
Instead, the HRDS nebulae probably span a wide range of evolutionary stages. Their small 
angular sizes result from their being located on average at large distances from the Sun. 



A. The HRDS Web Site 

We have constructed a web site to give others access to the HRDS datfl On this site 
one can download all the data from Tables [H [21 and [31 in addition to the spectra shown in 
Figure [21 the continuum scans shown in Figure [31 and the three-color MIR images shown in 
Figure [H We will continue to extend this site as more data are taken on the HRDS sources. 

The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. We would 
like to thank the anonymous referee whose careful reading has increased the clarity of this 
manuscript. The HRDS was partially supported by NSF award AST 0707853 to TMB. LDA 
was partially supported by SNF and by the NSF through GSSP awards 08-0030 and 09-005 
from the NRAO. LDA also acknowledges support from the ANR Agence Nationale award 
number ANR-08-BLAN-0241. This research has made use of NASA's Astrophysics Data 
System Bibliographic Services and the SIMBAD database operated at CDS, Strasbourg, 
France. 

Facility: Green Bank Telescope 



^http://go. nrao.edu/hrds 
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Table 1. GBT HRDS Source Catalog 



Name 


i 


6 


RA(J2000) 


Dec.(J2000) 


CH3OH'' 


CS 


Morphology'' Note'^ 




deg. 


deg. 


hh:mm:ss.s 


dd:mm:ss 


km 


km s"-*- 




GOOO.382+0.017 


0.382 


+0.017 


17:46:27.6 


-28:36:04 




84.4 


C 


GOOO.729-0.103 


0.729 


-0.103 


17:47:45.0 


-28:22:00 






I 


GOOO.838+0.189 


0.838 


+0.189 


17:46:52.1 


-28:07:20 


3.5 


5.4 


PB 


GOOl.324+0.104 


1.324 


+0.104 


17:48:20.2 


-27:45:02 






PB 


GOOl.330+0.088 


1.330 


+0.088 


17:48:24.9 


-27:45:12 






B 


GOOl. 330+0. 150 


1.330 


+0.150 


17:48:10.2 


-27:43:18 






IB {CN25) 


GOOl.488-0.199 


1.488 


-0.199 


17:49:53.6 


-27:45:54 






C 


G002.009-0.680 


2.009 


-0.680 


17:52:57.8 


-27:33:49 






B {CN28) 


G002.404+0.068 


2.404 


+0.068 


17:50:58.5 


-26:50:32 






C 


G002.418-0.611 


2.418 


-0.611 


17:53:37.8 


-27:10:34 






I 



"Meth anol maser velociti es are from: ICasw ell et alj l|l995'); 'EUingsen et al.' ('1996'); 'Gavlard & MacLeod' (fl993); 



MacLeod & Gavlard ( 1992'): iMentenI lll99l|) : Schutte et al, (1993); Slvsh ct al. (1999 ); Szvmczak ct al... (.2000,. 20021 ; 
van der Walt et al.l ( 119951) and lWalsh et al.i (Il997l . ll998h . 



'^Classification of source structure as seen in Spitzer GLIMPSE 8/im emission: 
B — Bubble: 8 /im emission surrounding 24 fira and radio continuum emission 

BB — Bipolar Bubble: two bubbles connected by a region of strong IR and radio continuum emission 

PB ~ Partial Bubble: similar to "B" but not complete 

IB - Irregular Bubble: similar to "B" but with less well-defined structure 

C - Compact: resolved 8 /ivci emission with no hole in the center 

PS - Point Source: unresolved 8 /tm emission 

I — Irregular: complex morphology not easily classified 

ND — Not Detected: source has no 8 fim emission 

We list in parenthesis the bubble identification name for sources in the IChurchwell et ah I hood boOTl ) IR bubble 
catalogs (see i|4.5l l. 

'^Notes on individual sources (see i) l4.4l l. PN? - Possible Planetary Nebula (see also ii l5.5ll . 



Note. — Table[T]is available in its entirety in the electronic edition of the Astrophysical Journal Supplement Series. 
A portion is shown here for guidance regarding its form and content. 
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Table 2. GBT HRDS Hydrogen Recombination Line Parameters 



Name 








I 




Tl 




AV 




Vlsr 




r.m.s. 








O 


C 




mK 


mK 


km s^-*- 


kms^^ 


km s^^ 


km s^-"- 


mK 


GOOO.382+0, 


.017a 


0, 


.382 


+0. 


,017 


137.6 


0.5 


23.0 


0.5 


25.7 


0.5 


2.4 


GOOO.382+0, 


.017b 





.382 


+0. 


.017 


104.5 


2.6 


32.4 


2.6 


41.4 


2.6 


2.4 


GOOO.729-0, 


.103a 


0, 


.729 


-0. 


.103 


136.2 


0.2 


22.6 


0.2 


105.3 


0.2 


2.9 


GOOO.729-0, 


.103b 


0, 


.729 


-0. 


.103 


114.4 


0.4 


33.1 


0.4 


83.2 


0.4 


2.9 


GOOO.838+0, 


.189 


0, 


.838 


+0. 


.189 


125.2 


0.1 


25.4 


0.1 


5.6 


0.1 


2.6 


GOOl.324+0, 


.104 


1, 


.324 


+0. 


.104 


104.7 


0.1 


27.0 


0.1 


12.7 


0.1 


2.4 


GOOl. 330+0, 


.088 


1, 


.330 


+0. 


.088 


100.2 


0.1 


27.4 


0.1 


13.2 


0.1 


2.7 


GOOl. 330+0, 


.150 


1, 


.330 


+0. 


.150 


30.6 


0.1 


21.7 


0.1 


20.1 


0.1 


2.4 


GOOl.488-0, 


.199 


1, 


.488 


-0, 


.199 


15.5 


0.3 


25.8 


0.3 


2.2 


0.3 


2.4 


G002.009-0, 


.680 


2, 


.009 


-0. 


.680 


17.3 


0.3 


30.2 


0.3 


18.2 


0.3 


2.1 


G002.404+0, 


.068 


2, 


.404 


+0. 


.068 


55.8 


0.1 


23.1 


0.1 


7.1 


0.1 


2.4 


G002.418-0, 


.611 


2, 


.418 


-0, 


.611 


8.8 


0.3 


18.6 


0.3 


7.8 


0.3 


2.0 



Note. — Table [2] is available in its entirety in the electronic edition of the Astrophysical Journal Supplement 
Series. A portion is shown here for guidance regarding its form and content. 
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Table 3. GBT HRDS Radio Continuum Parameters 



Name 


I 


b 






Ts 




ec 


CT0Q 






S 


aS 


Note'' 







O 


mK 


mK 


mK 


mK 


II 


II 


II 


II 


mJy 


mJy 




GOOO.382+0.017 


0.382 


+0.017 


2165 


14 


2239 


9 


105 


2 


88 


1 


1513 


53 


P 


GOOO.729-0.103 


0.729 


-0.103 


2335 


14 


2354 


7 


115 


1 


106 


1 


2131 


55 




GOOO.838+0.189 


0.838 


+0.189 


1025 


10 


1020 


8 


115 


2 


135 


2 


1176 


73 




GOOl.324+0.104 


1.324 


+0.104 


846 


5 


671 


8 


156 


2 


117 


2 


1035 


57 




GOOl. 330+0.088 


1.330 


+0.088 


869 


5 


704 


9 


161 


2 


105 


2 


991 


51 




GOOl. 330+0.150 


1.330 


+0.150 


137 


4 


145 


4 


128 


5 


99 


3 


134 


21 


P 


GOOl. 488-0.199 


1.488 


-0.199 


168 


5 


171 


4 


131 


5 


135 


6 


223 


37 


P 


G002.009-0.680 


2.009 


-0.680 


238 


4 


226 


5 


105 


3 


86 


3 


157 


19 


P 


G002.404+0.068 


2.404 


+0.068 


374 


8 


390 


19 


104 


3 


99 


7 


292 


67 


P 


G002. 418-0.611 


2.418 


-0.611 


80 


3 


78 


3 


104 


5 


101 


5 


61 


12 


P 



''Comments concerning continuum emission morphology and data quality. 
C - Complex: Source has complex continuum structure; this single Gaussian component model only crudely represents the 
true source characteristics (see ii l3.2ll . 

P - Peaked: Continuum peak lies within 10" of the nominal target position. These are the highest quality continuum data. 

Note. — Table [3] is available in its entirety in the electronic edition of the Astrophysical Journal Supplement Series. A 
portion is shown here for guidance regarding its form and content. 



-44- 



Table 4. HRDS Targets not Detected in RRL Emission 



Name 




Morphology^ 


Comments 


Object Type 


G003, 


,330-0, 


,870 


PS 


Compact IR and radio 


Possible PN 


G005, 


,756-0, 


,932 


B 


IR bubble with coincident radio 


Unknown 


G008, 


,362-0, 


,623 


PS 


Compact IR, coincident faint radio 


Likely PN 


G009, 


,741+0, 


,842 


PS 


Compact IR, coincident faint radio 


Likely PN 


GOll. 


,815+0, 


,842 


I 


Extended bright IR, coincident faint radio 


Likely Hll 


G012. 


,697+0, 


,683 


I 


Extended IR, coincident faint radio 


Likely Hll 


G013, 


,751-0, 


,360 


PS 


Compact IR and radio 


Possible PN 


G017, 


,223+0, 


,397 


PS 


Compact IR and radio 


Evolved star or PN 


G018, 


,446-0, 


,900 


PS 


Compact IR, faint radio 


Possible PN 


G019, 


,236+0, 


,495 


c 


Compact IR, relatively bright radio 


Possible PN 
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Extended IR, Faint radio emission 


Likely H II 
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B(N78) 


IR bubble, faint radio emission 
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Poor match between compact IR and diffuse radio 


Chance alignment 
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Poor match between compact IR and compact radio 


Likely AGN 
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,890 


ND 


Extended MIPSGAL, no GLIMPSE, compact radio 


Possible AGN 


G052, 
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Bright compact IR, compact radio 


Possible AGN 


G054, 


,096+0, 


,266 


I 


Pulsar Wind Nebula associated with SNR 54.1+0.3 


PWN 
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,161-0, 


,298 


C 


Extended GLIMPSE, faint radio 


Likely H II 


G064, 


,827+0, 


,494 


PS 


Poor match between compact IR and compact radio 


Likely AGN 
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,376+0, 


,280 


PS 


Compact IR and radio 


Possible PN 
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,386+0, 


,011 


ND 


Extended MIPSGAL, no GLIMPSE, faint radio 


Evolved star or PN 
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,081-0, 


,619 


PS 


Compact IR, faint radio 


Possible PN 



^Source structure classification as in Table [T] 
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Fig. 1. — Multi- wavelength infrared and radio images of four HRDS targets. From left to 
right are images made from Spitzer GLIMPSE 8.0 fim, Spitzer MIPSGAL 24 //m, and NRAO 
VGPS 21 cm continuum data. The HRDS target is at the center of each image. All panels 
are 5' on a side. The 82" GET resolution (FWHM beam size) is shown as a hatched circle. 
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Fig. 2. — Composite ( H n a ) RRL discovery spectra for the four Figure[T]H II regions. Shown 
is the average of the H 87 a through H 93 a RRLs, re-sampled to the resolution of the H 87 a 
spectrum and smoothed to 1.86 kms~^ resolution. A Gaussian fit to each hydrogen RRL 
component is superimposed. A vertical dashed line flags the nebula's hydrogen RRL LSR 
velocity, which is listed at the bottom of the flag. The expected locations of the helium and 
carbon RRLs are also flagged. The G031. 470— 0.344 spectrum clearly shows a helium line. 
For G030. 378+0. 106, hydrogen RRL emission is seen at two velocities. The 102.6 km s~^ 
component is probably associated with the W43 star forming complex (see §4.2p . 
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Fig. 3. — Radio continuum measurements of the four Figure[T]H II regions. Shown are the RA 
(top panels) and Dec (bottom panels) position scans centered at the nominal source position. 
Gaussian fits to the continuum emission that we associate with the nebula are superimposed. 
These typical data show that the continuum emission in the inner Galactic plane is often 
extremely complex. It is not always possible to distinguish between real continuum emission 
structure and sky noise. 
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Fig. 4. — Examples of HRDS target morphological classifications (see Tabled]). Each panel 
shows a three-color image constructed from Spitzer infrared data: MIPSGAL 24 /zm (red), 
GLIMPSE 8.0 /im (green) and GLIMPSE 3.6 /im (blue). The HRDS morphological classifi- 
cation is based on the target's appearance in the 8.0 fim image. Each image is 5' square and 
the 82" GBT beam is shown in the lower left corner. It is evident that the 8.0 /im emission 
always extends beyond that of the 24 /im emission. 
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Fig. 5. — RRL intensity as a function of angular distance from W43. Shown is the intensity 
of the 34 HRDS RRL components within 2° of the nominal center position of W43, {i, b) 
= (30?78, — 0?03) , whose velocities lie within the W43 velocity range, 85 to 105 kms~^ 
( iBalser et al.ll200l[ ). plotted as a function of the angular distance from the nominal center 
position. Error bars show the ±3 a uncertainties in Gaussian fits to the hydrogen RRL 
spectra. The decrease in RRL intensity with distance suggests that many of these multiple 
velocity components are associated with the W43 star forming complex. 
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Fig. 6. — Distribution of H II regions in Galactic longitude. Histograms show the distribution 
for the HRDS (dark gray) and the Known sample nebulae (dotted outline). 
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Fig. 7. — Distribution of H II regions in Galactic latitude. Histograms show the distribution 
for the HRDS (dark gray) and the known sample nebulae (dotted outline). Here we include 
only nebulae that are located within 1 ° of the Galactic plane. Gaussian fits to the two 
histograms are superimposed. The two samples have the same ~ 0?6 FWHM distribution 
in latitude. The HRDS source distribution is centered on the Galactic plane (— 0?01 is a 
displacement from b = 0° of less than 1/2 the GBT's HPBW), whereas the center of the 
known H II region sample distribution lies ~ 5' to the South. 
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Fig. 8. — Distribution of H II region radio recombina t ion lin e FWHM line widths. Histograms 
are shown for the HRDS (dark gray) and iLockman fll989l ) (dotted outline) H II region sam- 
ples. On average the HRDS nebulae have narrower lines than the Lockman sources. There 
are, moreover, many fewer HRDS nebulae with line widths greater than 35 kms~^. 
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Fig. 9. — Hll region RRL line intensity as a function of FWHM line width. Shown are the 
HRDS (left panel) and L89 samples (right panel). Vertical dotted lines mark line widths of 
35 kms~^. Low intensity HRDS sources share the same distribution of line widths as the 
entire sample. The L89 sources, however, show that systematically broader lines are found 
for low intensity sources. Vertical dashed lines mark line widths of 15 kms~^. Nebulae with 
line widths this narrow are "cool"; they cannot have electron temperatures greater than 
5,000 K. 
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Fig. 10. — Distribution of continuum angular size for HRDS nebulae. Histogram shows the 
geometric mean of the RA and Dec angular sizes derived from Gaussian fits to continuum 
scans (see Figure [3]) The vertical dotted line flags the 82 " HPBW size of the GET beam at 
X-band. The vertical dashed line flags the 116" size that would be measured by the GET 
for a source that just fills its beam. 
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Fig. 11. — Integrated source flux distributions for HRDS nebulae. The bottom panel shows 
the observed integrated continuum flux distribution. The top panel shows the integrated 
continuum flux distribution derived from the observed RRL line intensities and the contin- 
uum angular sizes, assuming a line-to-continuum ratio of 0.09. Solid lines are power law fits 
to each histogram. Dotted lines flag the point where the observed distribution deviates from 
the power law. Dashed lines show the power law fitted to the data in the other histogram. 
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Fig. 12. — Flux of H II regions ionized by single OB stars as a function of distance and spectral 
type. The gray curves show the estimated flux for Hll regions with electron temperatures in 
the range 5,000 to 10,000 K. The dashed line shows the > 70 mJy 9 GHz extrapolated flux 
limit of the HRDS target list, whereas the dash-dot line shows the > 180 mJy flux limit 
derived from the survey itself (see § 15. 6p . This flux limit implies that the HRDS is complete 
for Hll regions ionized by single 0-stars within 12kpc of the Sun (dotted line). 



